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ABSTRACT 


Various  meteorological  parameters  were  examined  to 
determine  their  predictive  value  for  forecasting  intensity 
changes  of  Western  North  Pacific  tropical  cyclones  using 
an  analog  forecast  scheme.   A  number  of  procedures  were 
devised  to  quantitatively  define  the  forecast  potentials 
of  these  parameters  in  a  systematic  method.   Based  on 
the  results  of  these  procedures,  predictors  for  the  24-hour, 
48-hour,  and  72-hour  forecast  periods  were  defined. 

To  evaluate  the  skill  of  the  predictors  and  the  pro- 
cedures used  to  define  them,  an  analog  intensification 
forecast  scheme  was  developed  and  applied  to  a  number  of 
test  cyclones.  The  results  of  these  forecasts  were  com- 
pared to  the  acceptability  criteria  established  by  the 
Joint  Typhoon  Warning  Center,  Guam  and  to  official  fore- 
casts issued  by  this  activity. 
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I.   INTRODUCTION 

A.   BACKGROUND 

Prior  to  the  mid-1960s,   studies  in  the  field  of  tropi- 
cal meteorology  progressed  at  a  very  slow  pace.   A  litera- 
ture search  through  various  meteorological  journals 
published  before  1965  reveals  a  strikingly  small  number  of 
articles  and  reports  pertaining  to  tropical  meteorology. 
However,  since  that  time,  the  amount  of  research  in  progress 
and  the  knowledge  of  the  tropical  ocean-air  regime  have  in- 
creased by  orders  of  magnitude.   Two  factors  have  been  pri- 
marily responsible  for  the  increased  efforts  and  results: 
the  meteorological  satellite  and  the  development  of  the 
efficient,  large  capacity,  second  generation  digital  com- 
puter.  Garstang  [1]  states  "...the  advent  of  the  meteoro- 
logical satellite,  and  the  power  of  the  computer  have 
collectively  led  to  great  effort  and  considerable  progress 
in  the  realm  of  tropical  meteorology...". 

The  one  phenomenon  of  tropical  meteorology  that  has  been 
the  object  of  more  research  than  any  other  is  the  tropical 
cyclone — the  hurricane  or  typhoon.   The  high  winds,  large 
storm  surges,  and  excessive  amounts  of  rainfall  associated 
with  these  systems  pose  a  highly  destructive  threat  to  many 
densely  populated  areas  of  the  world,  as  well  as  to  many 
high  density  shipping  routes  and  areas  of  extensive  naval 
operations.   Even  though  tropical  cyclones  lose  their 


organized  circulation  rapidly  as  they  encounter  large  land 
masses,  the  unstable  remnants  of  these  systems  can  cause 
wide-spread  damage  far  inland.   Their  far-ranging  destruc- 
tive nature  has  been  the  impetus  for  extensive  research  on 
the  natural  forces  which  initiate  and  maintain  tropical 
cyclones . 

Much  of  the  research  associated  with  tropical  cyclones 
has  been  directed  at  improving  the  accuracy  and  timeliness 
of  tropical  cyclone  forecasts.   The  major  factors  associated 
with  the  forecasting  are  the  initial  development,  movement 
and  intensification.   With  the  utilization  of  high  speed 
computers,  many  objective  movement  forecast  schemes  have 
been  developed  and  are  in  operational  use  as  forecast  aids. 
For  example,  there  are  currently  six  computerized  objective 
forecast  aids  for  typhoon  movement  in  use  at  Fleet  Weather 
Central/ Joint  Typhoon  Warning  Center  Guam  '  ( JTWC)  [ 2 ] .   None 
of  these  objective  schemes  are  considered  as  an  adequate 
forecast  technique,  but  are  used  as  aids  and  are  component 
of  the  final  product,  the  official  forecast.   However,  there 
are  currently  no  computerized  objective  schemes  in  opera- 
tional use  as  aids  to  forecasting  tropical  cyclone  intensi- 
fication.  The  bulk  of  all  forecast  research  efforts  has 
been  directed  at  movement,  for  if  the  movement  cannot  be 
predicted  adequately,  then  even  a  perfect  intensity  forecast 
may  have  limited  value.  The  overwhelming  emphasis  placed  on 
movement  rather  than  intensity  forecasts  is  documented  by 


the  fact  that  verification  of  the  latter  was  not  attempted 
at  JTWC  prior  to  1969.   But  with  the  increasing  accuracy  of 
movement  forecasts  that  has  been  experienced  in  recent 
years,  the  intensity  forecast  is  becoming  more  relevant  [3]. 

Due  to  the  lack  of  objective  aids  for  forecasting  tropi- 
cal cyclone  intensification,  the  forecaster  is  forced  to 
rely  on  the  limited  amount  of  current  and  recent  past  data 
at  hand  to  arrive  at  a  forecast.   At  JTWC  the  basic  forecast 
procedure  in  use  is  a  linear  extrapolation  of  past  intensity 
changes,  subjectively  modified  by  expected  conditions  along 
the  forecasted  track  and  by  climatology  [2].   These  proce- 
dures are  limited  in  that,  for  newly  developed  systems,  ex- 
trapolation is  especially  unreliable  (if  at  all  possible) 
and  rapid  deepening  cannot  be  forecast  to  an  adequate  degree 
of  accuracy.   An  objective  scheme  for  use  as  an  intensifica- 
tion forecast  aid  would  be  a  valuable  asset  to  the  tropical 
forecaster . 

One  of  the  objective  aids  for  typhoon  movement  employed 
at  JTWC  is  an  analog  scheme  called  TYFOON.   Forecasts  are 
provided  in  the  form  of  probability  ellipses  of  future  posi- 
tion,derived  from  the  movement  of  a  number  of  reports  of  past 
cyclones  that  have  been  screened,  selected  as  analogs,  and 
weighted  on  the  basis  of  their  similarity  to  the  subject 
current  cyclone.   The  degree  of  analog  similarity  is  a  func- 
tion of  time  of  occurrence  (month  and  day),  geographic  posi- 
tion, and  a  number  of  meteorological  parameters,  hereafter 
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called  analog  elements.   The  numerical  values  associated 
with  the  elements  of  a  particular  analog  are  compared  to  the 
respective  analog  elements  of  the  current  cyclone  to  deter- 
mine their  resemblance.   The  degree  of  resemblance  is  used 
to  assign  a  weight  to  the  respective  analog  which  will 
define  its  importance  to  the  final  results.   A  full  descrip- 
tion of  the  original  development  of  TYFOON  is  given  in  [4]. 
The  overall  results  of  TYFOON  have  been  the  best  of  any  ob- 
jective scheme  in  use  at  JTWC  [2],   The  lack  of  objective 
forecast  intensification  aids  and  the  success  of  the  TYFOON 
analog  scheme  for  movement  were  the  prime  stimulants  for 
initiating  this  study. 

B.   OBJECTIVES  OF  THE  STUDY 

The  development  of  the  TYFOON  forecasting  aid  necessi- 
tated the  selection  of  various  parameters  to  describe  the 
time,  position,  and  synoptic  features  of  current  and  past 
tropical  cyclones  in  the  western  North  Pacific  Ocean.   The 
historical  data  were  assembled  and  compiled  into  a  digital 
format  for  computer  utilization.   The  primary  objective  of 
this  study  was  to  conduct  a  pilot  study  of  the  analog  para- 
meters associated  with  TYFOON  to  determine  their  potential 
as  predictors  of  tropical-cyclone  intensification.   With  the 
assumption  that  some  positive  results  would  become  apparent 
from  the  primary  objective,  the  secondary  objective  was  to 
devise  a  forecast  scheme  that  would  test  the  predictive 
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value  of  the  analog  data  and  serve  as  a  first  attempt  at 
providing  the  forecaster  with  an  objective  analog  forecast- 
ing aid  to  predict  tropical-cyclone  intensification. 

C.  AREA  OF  STUDY 

The  area  of  study  was  restricted  to  the  region  for  which 
the  necessary  tailored  data  have  been  assimilated  and  digi- 
tized for  use  by  TYFOON .   This  area  consists  solely  of  the 
western  North  Pacific  Ocean.   Specifically,  it  is  that  region 
of  the  Pacific  north  of  the  equator  and  west  of  180  degrees. 
This  area  experiences  more  than  twice  the  number  of  tropical 
cyclone  occurrences  in  comparison  with  any  other  area  and, 
further,  is  the  only  region  where  the  systems  appear  in  all 
12  months  of  the  year  [5].   Consequently,  the  area  is  more 
suited  to  analog  forecasting  than  any  other  area  of  tropical 
cyclone  occurrence  due  to  the  availability  of  large  amounts 
of  tropical  cyclone  data. 

D.  DATA 

1 .   General 

The  data  available  for  the  study  were  dictated  by 
that  information  which  has  been  compiled  for  use  by  TYFOON. 
It  consists  of  various  parameters  of  two  major  categories: 
cyclone  identity   elements  and  analog  identity  elements. 
These  parameters  are  listed  in  Appendix  A.   The  data  have 
been  collected  and  digitized  for  each  six-hourly  report  of 
every  known  tropical  cyclone  (depression,  storm,  or  typhoon) 
since  1945,  and  at  the  time  of  this  study,  such  data  were 
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available  through  1969.   The  information  includes  best-track 
data,  a  result  of  post-season  analysis  of  operational  data. 
The  originally  observed  and  analyzed  data  have  been  re-ana- 
lyzed to  insure  consistency  in  development  and  movement  to 
a  degree  that  may  not  have  been  possible  at  the  time  of  col- 
lection . 

Only  limited  information  was  available  on  the  errors 
inherent  in  the  data.   Many  of  the  parameters  were  assigned 
numerical  values  as  a  result  of  subjective  analysis,  and 
therefore,  it  was  very  difficult  to  specify  the  expected 
errors  in  these  data.   It  is  believed  that  the  maximum  sus- 
tained surface  wind  has  an  error  of  10  to  15  percent  at  the 
time  of  observation  and  warning,  while  that  error  may  be 
reduced  to  10  percent  after  post  analysis  [6].   Estimates  of 
errors  associated  with  other  data  have  been  provided   and 
are  listed  in  Appendix  A. 

When  the  analog  elements  for  use  in  TYFOON  were 
selected,  only  cyclone  movement  was  being  considered. 
Therefore,  the  data  did  not  include  many  parameters  that  are 
closely  associated  with  typhoon  intensity  changes,  such  as 
thermal  structure  of  the  cyclone  position  of  the  monsoon 
trough,  or  sea-surface  temperature,  among  others.   However, 
there  were  some  parameters  such  as  the  past  12-hour  change 
in  sea-level  pressure,  the  past-movement  of  the  cyclone, 
and  the  minimum  700-mb  height  above  the  vortex,  which  are 


From  informal  discussion  with  CDR.  J.  D.  Jarrell. 
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known  to  have  some  relation  to  intensity  change,  at  least  in 
some  situations. 

2.   Selection  For  the  Study 

For  the  purpose  of  this  study,  several  test  cases 
were  selected,  and  for  each  test  case  a  number  of  analog 
candidates  were   defined  and  compared  to  the  test  case  with 
respect  to  intensification  characteristics  and  analog  element 
behavior.   A  test  case  was  defined  as  one  six-hourly  report 
of  a  tropical  cyclone.   An  analog  candidate  was  defined  as 
one  six-hourly  report  of  any  tropical  cyclone  (other  than 
the  one  from  which  the  test  case  was  selected)  which  by 
virtue  of  its  position  and  time  of  occurrence  (month  and  day) 
was  determined  to  be  similar  to  the  test  case. 

It  is  a  well  documented  fact  that  the  rate  of  occur- 
rence and  the  general  behavior  of  tropical  cyclones  vary 
from  season  to  season  and  month  to  month.   As  to  the  varia- 
tions in  intensification  behavior,  [7]  and  [8]  provide  sub- 
stantial evidence.   Since  this  study  was  aimed  at  the 
relationships  between  analog  elements  and  intensification, 
it  was  desirable  to  eliminate  other  variables  related  to 
intensification  (e.g.  month  of  occurrence)  without  severely 
restricting  the  number  of  possible  test  cyclones  and  analog 
candidates.   Therefore,  it  was  decided  to  select  for  the 
study  only  those  tropical  cyclones  occurring  in  one  season 
of  the  year.   It  was  also  considered  desirable  to  initiate 
the  research  using  somewhat  well-behaved  vorticies.   That 
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is,  a  season  in  which  there  occurs  relatively  little  rapid 
deepening  and  little  unsuspected  weakening  over  the  open 
ocean  was  considered  most  favorable  for  this  pilot  study. 
However,  it  was  decided  that  the  season  selected  should 
also  contain  a  relatively  large  number  of  occurrences,  for 
it  has  been  observed  that  better  results  are  obtained  from 

the  TYFOON  analog  scheme  when  larger  numbers  of  analogs  are 

2 
available.    Again,  it  was  felt  that  the  most  ideal,  yet  not 

overly  restrictive,  conditions  should  be  studied  initially. 

Studies  by  Brand  [7]  indicate  that  August  tropical 

cyclones  would  best  fulfill  the  above  requirements.   They 

exhibit  relatively  little  rapid  deepening,  while  also  showing 

little  weakening  over  low-latitude  open  ocean  areas.   August 

3 
also  has  the  second  highest  frequency  of  occurrence   of 

tropical  cyclones.   Therefore,  it  was  decided  to  select  at 
least  the  bulk  of  the  test  cases  from  this  month.   To  allow 
for  a  variable  time  span  from  which  analog  candidates  co  Ld 
be  selected  for  any  August  test  case,  it  was  necessary  t 
include  the  tropical  cyclones  of  July  and  September  in  t  ie 
data . 

Based  on  the  above  considerations,  the  analog  data 
for  the  months  of  July,  August  and  September  for  the  years 
1945  to  1969  were  obtained  (on  magnetic  tape)  from  the 
Environmental  Prediction  Research  Facility,  Monterey,  Calif- 
ornia . 


2 

From  informal  discussion  with  CDR.  J.  D.  Jarrell. 

3 

Number  of  six-hourly  reports. 
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II.   PROCEDURES 

A.  DEFINITION  OF  CYCLONE  INTENSITY 

In  a  study  of  tropical  cyclone  intensification,  the 
first  necessary  step  is  to  define  intensity.   The  variables 
most  commonly  associated  with  the  various  stages  of  develop- 
ment of  a  tropical  cyclone  from  depression  to  storm  to 
typhoon  or  hurricane  include  the  amount,  orientation,  and 
maximum  height  of  clouds,  the  minimum  sea-level  pressure, 
and  the  maximum  wind.   The  content  of  the  available  data 
eliminated  the  cloud  conditions  as  the  intensity  parameter 
for  this  study.   Of  the  remaining  two  parameters,  minimum 
sea-level  pressure  was  estimated  to  be  subject  to  smaller 
observational  errors  than  maximum  wind,  and  thus  would  be 
a  more  reliable  indicator  of  intensity.   However,  it  is  the 
maximum  wind  that  is  of  prime  interest  to  the  forecaster 
and  recipients  of  his  services;  therefore,  it  was  this  para- 
meter that  was  selected  to  describe  the  tropical  cyclone 
intensity . 

B.  SELECTION  OF  TEST  CYCLONES 

For  the  initial  study,  a  total  of  30  test  cases  were 
selected.   The  selection  was  not  random,  but  was  carried 
out  in  such  a  manner  to  insure  that  a  wide  range  of  initial 
intensities  and  positions  were  represented.   Widely  varying 
intensification  characteristics,  to  include  development  and 
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weakening,  were  also  exhibited  by  the  test  cases.   As  men- 
tioned before,  it  was  not  desirable  to  be  overly  restrictive 
with  the  data  as  this  would  possibly  inhibit  favorable  re- 
sults.  Based  on  the  findings  in  this  study,  restrictions 
could  be  applied,  if  necessary,  in  subsequent  studies.   In 
other  words,  if  one  adequate  objective  forecast  scheme  could 
not  be  derived  for  cyclones  of  widely  ranging  initial  in- 
tensities and  positions,  then  the  input  data  could  be  strati- 
fied with  respect  to  these  variables  and  different  forecast 
schemes  attempted  for  the  various  classes  of  cyclones. 

Other  factors  were  considered  in  the  test  cyclone  selec- 
tion.  Since  one  avenue  of  research  utilized  historical  data 
of  the  test  cases,  a  restriction  was  placed  on  the  minimum 
amount  of  history  possessed  by  a  number  of  the  cases.   There 
were,  however,  some  test  cases  with  no  recorded  historical 
data.   Recurvature  was  also  considered  and  a  number  of  the 
test  cyclones  were  selected  on  the  basis  of  whether  or  not 
recurvature  took  place  within  the  forecast  periods  or  in  the 
recent  past . 

As  stated  before,  the  bulk  of  the  test  cases  were 
selected  from  the  month  of  August.   However,  to  insure  that 
a  sufficient  number  of  appropriate  test  cases  from  recent 
years  collectively  represented  a  wide  range  of  initial  in- 
tensities and  positions,  and  to  avoid  an  excessive  number 
of  related  test  cases  (from  the  same  tropical  cyclone)  it 
was  found  necessary  to  select  a  few  test  cyclones  from  the 
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month  of  September.   The  30  test  cases  initially  selected 
for  the  study  are  listed  in  Appendix  B. 

The  total  number  of  six-hourly  reports  included  in  the 
three-month  period  of  data  obtained  for  the  study  amounted 
to  over  7500.   It  is  realized  that  only  30  test  cases  may 
not  constitute  a  representative  sample  of  all  analog  para- 
meters involved.   However,  the  full  range  of  possible 
initial  intensities  and  a  highly  variable  distribution  of 
initial  positions  were  represented.   It  was  believed  that 
if  useful  persistent  relationships  (between  the  analog 
elements  and  intensification)  not  only  exist,  but  also  re- 
main valid  for  highly  variable  initial  conditions,  then  at 
least  some  indications  should  be  observed  from  as  few  as 
30  test  cases.   Again,  it  is  emphasized  that  the  prime 
purpose  of  this  study  was  to  determine  if  some  analog 
elements  possess  any  potential  as  predictors  of  intensifi- 
cation, not  to  quantitatively  describe  the  various  degrees 
of  this  potential. 

C.   BASIC  APPROACH  TO  THE  MAIN  OBJECTIVE 

The  main  objective  of  this  study  required  the  develop- 
ment of  a  procedure  that  would  examine  the  analog  elements 
and  identify  any  behavior  which  appeared  to  be  related  to 
intensification.   Several  trial  and  error  forecast  schemes 
could  have  been  attempted  using  various  parameters  with 
associated  weight  factors  being  assigned  intuitively.   It 
was  decided,  however,  that  a  systematic  approach  would  be 
more  likely  to  bear  fruitful  results. 
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In  an  analog  intensification  forecast  scheme,  a  compari- 
son would  be.  made  of  known  elements  (predictors)  of  the 
current  cyclone  with  those  of  past  cyclones  that  occurred 
within  specified  space  and  time  limits,  as  well  as  exhibiting 
an  intensity  within  some  predetermined  range.   Based  on  a 
comparison  of  the  analog  elements,  a  number  of  analog  candi- 
dates that  most  closely  resemble  the  current  cyclone  would 
be  determined.   The  changes  in  intensity  exhibited  by  these 
most  similar  analogs  would  be  formulated  into  a  forecast  for 
the  current  cyclone. 

To  investigate  the  forecast  potential  of  the  analog  ele- 
ments, it  was  decided  to  use  the  inverse  of  an  analog  fore- 
cast procedure.   For  a  given  test  case  of  known  future 
intensification,  all  six-hourly  reports  of  past  cyclones 
were  screened  to  determine  those  that  occurred  within  speci- 
fied time,  space,  and  initial  intensity  limits.   Then, 
rather  than  classify  the  analogs  on  the  basis  of  their 
element  values  (as  would  be  done  in  a  forecast  scheme)  the 
future  intensification  of  each  selected  analog  candidate  was 
determined,  and  they  ware  classified  according  to  their 
deviations  in  future  intensification  from  that  of  the  cur- 
rent cyclone.   A  class  was  defined  as  a  collection  of  analogs 
that  exhibited  similar  deviations  in  intensification  from 
that  of  the  test  case  for  a  given  forecast  period.   Each 
class  of  analogs  represented  an  intensification  range,  and 
the  similarity  of  that  range  to  the  test  cyclone  intensifi- 
cation varied  from  class  to  class.   The  analog  elements 
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associated  with  the  candidates  in  each  class  were  ranked  to 
determine  their  average  degree  of  matching  with  respect  to 
the  test  case  elements.   This  basic  approach  will  become 
better  understood  as  the  various  procedures  are  discussed. 

The  hypothesis  around  which  the  basic  approach  to  the 
study  was  designed  specifies  that  if  any  of  the  analog  ele- 
ments were  related  to  intensification,  then  the  degree  of 
matching  of  these  elements,  with  respect  to  the  current 
cyclone  values,  should  decrease  as  the  deviation  in  inten- 
sification of  the  associated  analog  candidates,  from  that 
of  the  test  case,  increases.   If  the  deviations  in  inten- 
sification are  thought  of  as  forecast  errors,  then  the 
hypothesis  implies  that  higher  degrees  of  matching  are 
associated  with  smaller  forecast  errors  for  potential  pre- 
dictors of  intensification.   If  any  of  the  elements 
exhibited  this  trend  with  some  degree  of  consistency,  then 
a  class  of  analogs  with  intensification  similar  to  the  test 
case  could  be  determined,  and  this  intensification  used  as 
a  forecast  aid.   This  basic  procedure  was  considered  to 
be  the  most  efficient  and  direct  approach  to  the  main  objec- 
tive . 

The  main  problems  encountered  in  this  basic  approach 
were  deciding  on  what  criteria  to  use  for  the  analog  classi- 
fication and  what  technique  to  use  for  the  degree  of  matching 
computation.   The  following  sections  describe  the  various 
specific  procedures  that  were  developed  in  an  attempt  to 
accomplish  the  primary  and  secondary  objectives.   All 
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procedures  were  programmed  in  FORTRAN  IV  for  use  with  the 
IBM  360  computer  at  the  Naval  Postgraduate  School,  Monterey, 
California . 

D.   FIRST  APPROACH  TO  THE  OBJECTIVE 

The  first  procedure  developed  for  the  study  was  an  in- 
verse analog  program,  as  described  in  the  previous  section. 
It  consisted  of  two  basic  operations:  (1)  test  case  defini- 
tion and  analog  selection;  (2)  analog  classification  and 
parameter  rank  computation. 

In  the  first  operation  a  test  cyclone  was  defined  by 
means  of  a  test  case  identification  number,  by  its  associ- 
ated latitude,  longitude,  and  month,  day  and  time  of  occur- 
rence, and  by  its  associated  analog  elements.   The  future 
intensifications  for  12,  24,  48,  and  72  hours,  exhibited 
by  the  test  case,  were  also  defined.   For  test  cases  not 
exhibiting  a  particular  future  intensification,  a  code  of 
999  was  assigned.   The  intensification  for  each  forecast 
period  was  defined  as  the  difference  between  the  intensit 
at  the  end  of  the  respective  period  and  the  initial 
intensity . 

The  next  step  was  to  define  the  time,  space,  and  initial 
intensity  limits  to  be  utilized  in  the  analog  selection. 
Time  and  space  screening  criteria  of  plus  or  minus  15  days 
from  the  day  of  the  test  cyclone,  and  plus  or  minus  five 
degrees  latitude  and  longitude  from  the  initial  position  of 
the  test  case  were  selected  for  this  procedure.   In  general, 
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vorticies  of  extremely  different  intensities  exhibit  unlike 
intensification  characteristics  over  an  extended  period  of 
time.   For  example,  a  young  developing  tropical  depression 
would  be  unlikely  to  exhibit  intensification  characteristics 
similar  to  that  of  a  fully-developed  super  typhoon.   Under 
conditions  favorable  for  development,  the  young  cyclone 
could  continue  to  intensify  for  a  considerable  length  of 
time,  whereas  the  super  typhoon,  being  near  or  at  its  point 
of  maximum  intensity,  would  not  be  expected  to  exhibit  this 
behavior.   In  an  analog  forecast  scheme,  an  initial  intensity 
screen  would  be  utilized  to  eliminate  unnecessary  and  un- 
desirable noise  in  the  selected  analogs.   Since  the  analog 
selection  phase  of  this  procedure  was  similar  to  that  which 
would  be  utilized  in  an  analog  forecast  scheme,  the  initial 
intensity  screen  was  applied.   Limits  of  plus  or  minus  20 
knots  from  that  of  the  test  case  were  defined  for  the  analog 
selection.   With  the  test  cyclone  and  its  associated  screen- 
ing criteria  established,  the  procedure  continued  with  the 
analog  selection. 

Each  six-hourly  report  of  all  past  tropical  cyclones 
(contained  in  the  data) ,  excluding  the  cyclone  with  which 
the  test  case  was  associated,  was  screened  with  respect  to 
its  time  of  occurrence  (month  and  day),  position,  and  initial 
intensity.   All  those  reports  satisfying  the  previously 
defined  screening  criteria  were  designated  analog  candidates. 
For  each  analog  selected,  the  associated  12-,  24-,  48-,  and 
72-hour  future  intensifications  were  determined.   With  the 
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analogs  selected  and  their  intensifications  computed,  the 
first  operation  in  this  procedure  was  completed. 

The  analog  classification  and  parameter  rank  computation 
proved  to  be  the  most  difficult  aspect  of  the  study.   No 
previous  guidelines  were  available  to  indicate  what  would  be 
the  best  approach.   The  various  methods  of  attack  were  some- 
what trial  and  error,  but  all  were  based  on  the  hypothesis 
stated  earlier,  and  all  had  one  concept  in  common,  namely 
the  analogs  were  classified  according  to  their  deviations 
in  intensification  from  that  of  the  test  case  for  each  fore- 
cast interval . 

In  the  first  approach  to  the  analog  classification, 

positive  and  negative  deviations  were  considered.   For  each 

forecast  period  the  following  analog  classification  scheme 

was  implemented: 

Class  1  Deviation  =    0->  +9  KTS 

Class  2  Deviation  =    0+  -9  KTS 

Class  3  Deviation  =   10+  19  KTS 

Class  A  Deviation  =  -10+-19  KTS 

•  •  • 

•  •  • 
Class  13     Deviation  =   70->  79  KTS 
Class  14     Deviation  =  -7CK-79  KTS 

With  all  (or  most)  of  the  analog  candidates  classified,  the 
next  step  was  to  compute  analog  element  ranks.   Analogs 
with  deviations  in  intensification  not  defined  above  were 
eliminated  from  further  investigation. 

For  each  analog  element  associated  with  an  analog  can- 
didate in  the  above  classes,  a  rank  was  computed  to  deter- 
mine the  degree  of  matching  with  the  respective  test  cyclone 
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element.   In  this  approach  the  ranges  of  the  analog  elements 


over  the  entire  data  file  were  utilized.    The  rank  of  the 
N    element  of  each  classified  analog  candidate  was  computed 
using  the  following  formula: 


R(N)  = 

PA(N)  > 
PT(N)  ■ 
r(N) 
R(N) 


(1.0  -  r*v";  ~  rx^;  )  100.0 
r  (N) 

N  .  element  value  of  the  analog 

t  h 
N    element  value  of  test  cyclone 

range  of  N*-"  element 

rank  of  Nfc"  element  associated  with  the 

particular  analog 


(1) 


The  ratio 


PA(N)  -  PT(N) 
r(N) 
th 


indicates  how  different  the  two 
values  of  the  N"   element  are  as  compared  to  the  maximum 
possible  difference.   Subtracting  this  quantity  from  1.0 
results  in  a  higher  rank  for  a  closer  match.   For  an  identi- 
cal match,  PA(N)  =  PT(N)  and  R(N)  =  100.0;  a  rank  of  zero 
would  indicate  the  two  values,  PA(N)  and  PT(N),  were  on 
opposite  extremes  of  the  range. 

Once  the  elements  were  ranked  for  all  classified  analo  s, 
the  average  rank  of  each  element  was  computed  for  each  cl<-  ;s, 
thus  arriving  at  an  average  degree  of  matching.   The  standard 
deviation  of  the  individual  ranks  in  each  class  was  also 
calculated  for  each  analog  element.   Also  for  each  element, 
a  graph  of  average  rank  versus  class  deviation  in  intensifi- 
cation was  plotted.   Recalling  that  these  deviations  could 


Maximum  and  minimum  values  of  the  analog  elements  and  the 
ranges  were  computed  in  a  separate  program  and  specified 
in  this  procedure.   For  direction  of  movement  elements, 
the  range  was  specified  as  180  degrees. 
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be  considered  forecast  errors,  these  graphs  were  essentially 
the  average  degree  of  matching  versus  forecast  error. 

This  procedure  was  applied  to  the  first  15  test  cyclones 
listed  in  Appendix  B,  and  the  results  were  analyzed.   Graphs 
revealing  a  decreasing  average  degree  of  matching  with  in- 
creasing forecast  error  were  not  observed  with  any  persis- 
tence for  any  of  the  analog  elements.   In  several  cases,  an 
increasing  average  degree  of  matching  with  increasing  fore- 
cast error  resulted  with  the  past  direction  and  speed  of 
movement  elements.   No  physical  interpretation  could  be 
given  for  this  behavior. 

One  obvious  fault  in  the  procedure  was  noted.   For  the 
48-  and  72-hour  forecast  periods,  and  to  a  lesser  extent  the 
24-hour  period,  the  selected  analog.s  in  any  one  class  ex- 
hibited markedly  different  intermediate  intensification 
characteristics.   In  a  48-hour  period,  for  example,  one  ana- 
log may  have  intensified  30  knots  in  a  gradual  but  steady 
manner,  while  another  analog  may  have  intensified  50  knots 
in  the  first  24  hours,  remained  steady  for  12  hours  and  then 
weakened  by  20  knots  in  the  next  12  hours.   These  two  hypo- 
thetical cases  exhibit  markedly  different  behavior,  but 
since  each  experiences  a  net  48-hour  intensification  of  plus 
30  knots  they  would  be  classified  together  under  the  current 
procedures.   Since  it  was  the  intention  to  classify  only 
similar  analogs  together  in  order  to  study  the  variable  be- 
havior of  the  analog  elements,  modified  definitions  of 
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intensification  were  implemented  for  the  24-,  48-,  and  72- 
hour  forecast  periods,  and  are  shown  below.   The  letter  I 
represents  intensity  and  the  subscripts  indicate  the  number 
of  hours  from  the  time  of  occurrence  of  the  test  cyclone 
or  analog . 

A  I24  =  T"24  "  I12 


A  X48  =  X48  "  X24  (2) 


A  X72  =  hi    ~   X48 


For  purposes  of  analog  classification  for  24-hour  intensi- 
fication, the  deviation  in  intensification  from  that  of  the 
test  case  between  time  =  12  and  time  <=  24  had  to  fall  within 
a  set  of  class  limits.   For  48-hour  intensification,  the 
deviations  for  each  12-hour  period  involved  (24  to  36  and  36 
to  48  hours)  had  to  fall  within  one  particular  set  of  class 
limits.   Similar  restrictions  were  placed  on  the  analog 
classification  for  the  72-hour  intensification.   Classifica- 
tion for  the  12-hour  period  remained  unchanged.   Compared 
to  the  original  definitions  of  intensification,  these  defini- 
tions would  require  more  time-restricted  relationships  to 
exist  between  the  analog  elements  and  intensification  in 
order  for  the  desired  trend  of  decreasing  degree  of  element 
matching  with  increasing  forecast  error  to  appear  with  any 
persistence.   As  an  example,  for  an  element  to  exhibit  a 
strong  relationship  to  48-hour  intensification,  the  relation- 
ship would  have  to  be  valid  during  more  precisely  defined 
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time  limits  (24  to  48  hours).   Under  the  prior  definitions 
of  intensification,  an  indicated  relationship  between  the 
behavior  of  an  element  and  the  48-hour  intensification 
could  have  been  actually  very  strong  for  the  first  24  hours, 
but  relatively  weak  for  the  second  24  hours.   However , these 
restrictions  were  considered  more  favorable  than  to  contend 
with  the  identical  classification  of  analogs  exhibiting 
highly  variable  intensifications. 

The  first  15  test  cases  were  rerun  using  the  modified 
definitions  of  intensification,  and  the  results  evaluated 
again.   Little  if  any  change  from  the  previous  results  were 
observed.   In  most  cases  the  average  ranks  of  all  the 
elements  were  very  high,  generally  greater  than  85.   This 
was  due  to  the  use  of  the  element  ranges  in  the  equation 
for  the  ranks.   The  computed  ranges  were  usually  the  result 
of  extreme  and  uncommon  element  values.   For  example,  the 
range  of  sea-level  pressure  extended  from  875  mb  to  1012  r  >. 
To  compensate  for  this,  the  standard  deviations  of  the 
analog  element  distributions  over  the  entire  data  file  we  e 
computed  and  substituted  for  the  ranges.   With  the  use  of 
the  standard  deviations,  negative  ranks  would  result  when 
the  difference  in  values  of  an  element  was  greater  than  the 
standard  deviation  of  the  element  distribution.   A  closer 
match  of  an  analog  element  to  its  test  cyclone  value  would 
result  in  a  more  positive  rank. 

A  third  run  of  the  first  15  test  cases  revealed  no 
positive  results.   Plots  of  the  average  rank  versus  forecast 
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error  revealed  a  much  more  variable  pattern  than  previously 
observed,  but  no  consistent  decreasing  degree  of  matching 
with  increasing  forecast  error  resulted  for  any  of  the  analog 
elements.   Many  cases  revealed  extreme  oscillations  of  the 
average  ranks  of  some  elements  from  one  class  of  analogs  to 
the  next . 

Only  one  relationship  persisted  throughout  the  entire 
study  to  this  point.   In  almost  every  instance  for  any  analog 
element  in  any  analog  class,  a  larger  average  rank  was 
associated  with  a  smaller  standard  deviation  of  the  element 
ranks  involved.   Only  a  small  increase  in  an  average  rank 
was  generally  accompanied  by  a  small  decrease  in  the  standard 
deviation.   Therefore,  lower  average  ranks  were  associated 
with  relatively  more  variable  individual  ranks,  but  higher 
average  ranks  were  associated  with  less  variable  individual 
ranks . 

The  fact  remained,  however,  that  both  high  and  low 
average  ranks  were  observed  in  all  classes  of  analogs  for 
all  15  test  cases.   It  was  felt  that  if  some  persistent  and 
useful  trend  were  to  develop,  at  least  some  weak  indications 
would  have  resulted  by  this  point  in  the  research.   In  the 
absence  of  any  such  indications,  this  method  of  analog 
classification  and  element  ranking  was  abandoned.   The  possi- 
bility existed  that  the  ranking  procedure  utilized  was  too 
sensitive  for  the  accuracy  of  the  data  involved  and  for  the 
reliability  and  persistence  of  some  relatively  weak,  but 


28 


somewhat  useful  relationships  which  might  exist.   It  was 
decided  to  attempt  a  less  restrictive  classification  pro- 
cedure and  to  adopt  a  less  refined  ranking  technique. 

E.   SECOND  APPROACH  TO  THE  OBJECTIVE 

The  two  basic  operations  of  the  first  approach  were 
also  utilized  in  this  procedure.   The  test  cyclone  defini- 
tion and  analog  selection  phase  was  modified  in  two  ways. 
First,  the  initial  intensity  limits  for  the  analog  screen- 
ing were  expanded  to  plus  or  minus  40  knots  from  the 
initial  intensity  of  the  test  case.   It  was  felt  that  a 
greater  number  of  analogs  would  generate  better  results,  as 
has  been  observed  with  the  TYFOON  analog  scheme.   Due  to 
possible  diurnal  effects  on  the  results,  the  12-hour  fore- 
cast was  excluded  and  the  24-hour  intensification  re- 
defined as 


A  X24  "  *24  "  Z0 


(3) 


The  analog  classification  and  element  rank  computation 
were  completely  re-developed  for  this  procedure.   For  analog 
classification,  the  deviations  in  intensification  during 
the  24-hour  periods  of  initial  time  to  24  hours,  24  to  43 
hours,  and  48  to  72  hours  were  considered.   The  deviations 
within  the  individual  12-hour  periods  contained  in  these  24- 
hour  periods  were  not  considered  in  order  to  make  the  classi- 
fication less  restrictive  and  to  eliminate  possible  diurnal 
effects.   For  each  forecast  period  (24-,  48-,  and  72-hours), 
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the  analogs  that  intensified  within  plus  or  minus  10  knots 
of  the  test  case  per  24-hour  period  were  identified  and 
classified  together. 

The  next  step  was  to  assign  a  rank  to  each  element  of 
the  classified  analogs.   Rather  than  define  a  degree  of 
matching  for  each  individual  element,  a  technique  was  de- 
veloped that  would  indicate  a  match  or  non-match  result 
for  each  element  and  all  matches  would  be  given  equal 
weight.   It  was  necessary  to  establish  for  each  analog 
element  an  interval,  centered  on  the  test  case  value,  within 
which  an  analog  candidate  value  would  have  to  occur  to  be 
considered  a  match.   If  a  match  occurred,  the  element  of 
the  analog  was  assigned  a  rank  of  1.0.   If  a  match  failed 
to  occur,  the  element  was  assigned  a  rank  of  0.0. 

Due  to  the  variable  ranges  of  possible  values  among 
the  analog  elements  each  established  interval  of  acceptance 
had  to  be  unique  to  a  particular  element.   Since  the 
standard  deviation  of  a  data  population  indicates  the  va  ia- 
bility  of  the  quantity,  the  interval  of  acceptance  was  ir  ide 
proportional  to  this  measure  of  dispersion.   An  interval 
of  plus  or  minus  one  fourth  of  a  standard  deviation  centered 
on  the  test  case  value  was  selected  as  the  interval  of 
acceptance  for  each  analog  element.   All  analog  element 
values  occurring  within  this  range  were  considered  a  match; 
all  others  were  rejected.   The  acceptance  intervals  for  the 
analog  elements  are  given  in  Appendix  C. 
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For  each  element  the  number  of  occurrences  of  missing 
data  associated  with  the  classified  analogs  was  determined. 
Subtracting  these  numbers  from  the  total  number  of  analogs 
in  the  class  resulted  in  the  number  of  possible  matches 
for  each  element.   By  summing  up  the  ranks  associated  with 
a  particular  element  over  all  the  classified  analogs,  the 
number  of  matches  was  determined.   Dividing  this  by  the 
number  of  possible  matches  and  multiplying  by  100  resulted 
in  the  percent  number  of  matches  for  the  element.   This 
quantity  was  now  referred  to  as  the  degree  of  matching. 
This  procedure  was  applied  and  a  degree  of  matching  estab- 
lished for  each  of  the  18  elements. 

The  classification  procedure  was  now  repeated  for  all 
the  analogs  that  exhibited  an  intensification  within  plus 
or  minus  20  knots  of  the  test  cyclone  per  24-hour  period, 
and  for  each  forecast  period.   The  degree  of  matching  of 
each  element  was  computed  exactly  as  before.  The  procedures 
were  again  repeated  for  deviations  within  plus  or  minus  30 
knots  and  finally  for  deviations  within  plus  or  minus  40 
knots  per  24-hour  period,  for  each  forecast  period. 

For  each  element,  the  degree  of  matching  versus  the 
maximum  allowable  class  deviation  in  intensification  (maxi- 
mum allowable  forecast  error)  was  plotted.   As  in  the  first 
approach  a  decreasing  degree  of  matching  with  increasing 
forecast  error  was  considered  to  be  indicative  of  a  relation- 
ship between  element  behavior  and  intensification. 

As  the  maximum  allowable  forecast  error  increased  from 
10  to  40  knots  the  number  of  analogs  and  therefore  the  number 
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of  possible  matches  increased.  But  for  a  potential  predic- 
tor, any  increased  number  of  actual  matches  was  expected  to 
result  in  a  lower  degree  of  matching. 

The  procedure  was  applied  to  all  30  test  cases  listed 
in  Appendix  B,  and  the  results  analyzed.   The  desired  trend 
of  decreasing  degree  of  matching  with  increasing  forecast 
error  was  exhibited  to  at  least  some  extent  by  all  analog 
elements  but  much  more  frequently  by  some  elements. 

The  next  step  was  to  determine  the  relative  signifi- 
cance of  each  analog  element  as  a  potential  predictor  of 
intensification.   For  each  forecast  period  and  for  each  of 
the  analog  elements,  the  number  of  occurrences  of  decreasing 
degree  of  match  with  increasing  forecast  error  was  noted  and 
the  percent  of  this  number  to  the  total  possible  number  of 
occurrences  was  computed.  These  percentages  were  referred 
to  as  indicated  forecast  potentials.   The  total  possible 
number  of  occurrences  varied  with  individual  test  cases. 
Some  cases  contained  missing  data  which  ruled  out  possible 
matches  for  some  elements  and  in  some  cases  a  future  48- 
or  72-hour  verifying  intensity  did  not  exist  which  prohibited 
analog  classification.   The  results  of  these  computations 
are  presented  in  Appendix  D. 

As  a  first  attempt,  it  was  decided  to  choose  for  each 
forecast  period  the  analog  elements  associated  with  the 
five  largest  indicated  forecast  potentials  as  predictors 
of  intensification  for  the  respective  period. 
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Before  developing  an  experimental  forecast  scheme,  one 
further  step  was  conducted  to  provide  a  more  rigorous  study 
of  the  analog  elements. 

F.   THIRD  APPROACH  TO  THE  OBJECTIVE 

In  order  to  make  the  study  as  thorough  as  time  would 
allow,  various  parameters  were  combined  to  derive  other  ana- 
log elements  to  be  examined  for  forecast  potential.   Also, 
since  it  was  the  change  in  intensity  that  was  to  be  forecast, 
the  recent  time  rates  of  change  of  some  original  and  derived 
elements  were  investigated  for  forecast  potential.   Only  the 
past  24-  and  48-hour  changes  of  various  parameters  were 
examined,  as  utilizing  any  greater  amount  of  history  might 
lead  to  limitations  in  a  forecast  scheme  applied  to  a  re- 
cently detected  cyclone.   The  derived  elements  are  listed  in 
Appendix  E.   All  but  three  were  easily  derivable  from  the 
original  analog  and  identify  elements  of  TYFOON. 

The  three  new  elements  which  required  some  external  in 
formation  were  the  1000/700-mb  thickness  elements.   These 
elements  were  included  as  indicators  of  the  thermal  charac- 
teristics of  the  lower  domain  of  the  cyclones.   The  minimum 
700-mb  height  above  the  cyclones  and  the  minimum  sea-level 
pressure  (original  analog  elements)  made  possible  an  approxi- 
mation of  the  1000/700-mb  thickness. 

The  missing  term  for  the  thickness  elements  was  the 
height  of  the  1000-mb  surface.   It  is  a  widely-held  belief 
that  a  sea-surface  temperature  of  at  least  26-27C  is  required 
for  tropical  cyclone  development  [5].   The  air  temperature 
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near  the  air-sea  interface  must  be  similar  to  the  sea-sur- 
face temperature.   A  mean  virtual  temperature  of  30C  for 
the  layer  between  the  surface  and  1000  mb  was  considered  a 
reasonable  average  approximation.   With  this  approximation, 
and  the  use  of  the  hydrostatic  equation  and  the  equation  of 
state,  the  height  of  the  1000-mb  surface  could  be  estimated 
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The  1000-mb  height  thus  obtained  was  subtracted  from  the 
minimum  700-mb  height  above  the  vortex  to  determine  the 
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thickness.   In  this  manner,  the  current,  past  24-hour  and 
past  48-hour  thicknesses  were  computed  and  the  appropriate 
time  rates  of  change  calculated. 

The  investigative  procedures  conducted  on  these  derived 
elements  were  essentially  identical  to  those  applied  to  the 
original  elements  and  described  in  the  previous  section. 
Some  program  modifications  were  necessary,  however. 

In  the  test  cyclone  definition  and  analog  selection 
phase,  the  necessary  past  24-  and  48-hour  data  were  defined 
and  a  routine  was  added  to  compute  the  time  rates  of  change. 
A  history  indicator  was  also  defined  for  each  test  cyclone 
to  indicate  the  extent  to  which  history  elements  of  the 
selected  analogs  should  be  calculated.   This  phase  was  also 
modified  to  compute  the  necessary  and  available  derived 
elements  for  each  selected  analog. 

Only  one  modification  was  necessary  for  the  analog 
classification  and  element  rank  computation  phase.   It  was 
necessary  to  redefine  the  matching  criteria  (acceptance 
intervals)  for  the  new  elements.   To  compute  the  standard 
deviations  of  all  the  elements  would  have  required  more  com- 
puter time  than  was  practically  available.   Acceptance  inter- 
vals were  defined  subjectively  based  on  the  estimated 
variability  of  each  element  and  on  the  estimated  probable 
errors  in  the  values.   As  before  the  intervals  were  centered 
on  the  test  cyclone  element  values,  and  for  an  analog 
element  to  be  accepted  as  a  match,  its  value  had  to  occur  in 
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the  appropriate  interval.   The  acceptance  intervals  for  the 
derived  elements  are  listed  in  Appendix  E. 

The  procedures  were  carried  out  on  the  30  test  cases 
listed  in  Appendix  B  and  an  analysis  of  the  results  was  made 
as  described  in  the  previous  section.   For  each  forecast 
period,  the  indicated  forecast  potential  was  computed  for 
each  element,  and  the  elements  associated  with  the  five  lar- 
gest indicated  forecast  potentials  were  selected  as  predic- 
tors of  intensity.   The  indicated  forecast  potentials  of 
the  derived  analog  elements  are  listed  in  Appendix  F. 

With  the  intensity  predictors  defined,  it  was  now  neces- 
sary to  establish  their  relative  importance  for  use  in  a 
forecast  scheme. 

G.   DERIVATION  OF  WEIGHT  FACTORS  FOR  SELECTED  PREDICTORS 

The  results  of  the  study  to  this  point  revealed  that  the 
potential  among  the  original  and  derived  analog  elements  as 
predictors  of  tropical  cyclone  intensification  varied 
greatly.   For  the  elements  that  were  selected  as  predict  >rs, 
in  a  forecast  scheme,  more  emphasis  would  be  placed  on  those 
that  exhibited  more  persistently  the  trend  of  decreasing  de- 
gree of  matching  with  increasing  forecast  error.   The  degree 
of  emphasis  placed  on  each  predictor  would  be  indicated  by 
a  weight  factor.   Out  of  ten  predictors,  for  example,  it 
would  not  suffice  to  assign  the  best  one  a  weight  of  ten  and 
the  worst  one  a  weight  of  one.   This  would  imply  that  the 
best  predictor  possesses  ten  times  the  forecast  potential  of 
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the  worst  one,  whereas  it  may  only  possess  twice  as  much. 
Therefore  it  was  necessary  to  derive  weight  factors  as  a 
function  of  the  indicated  forecast  potentials  resulting 
from  the  earlier  investigations. 

To  assign  weight  factors  relative  to  the  indicated 
results,  the  following  procedures  were  applied.   First  of 
all,  a  small  number  of  selected  predictors  were  eliminated 
due  to  their  close  relationships  to  other  predictors  with 
equal  or  higher  indicated  forecast  potentials.   This  avoided 
a  concentrated  emphasis  on  the  related  predictors,  the  cumu- 
lative effect  of  which,  might  often  lead  to  poor  results. 
For  each  forecast  period  the  selected  predictor  (original 
or  derived)  with  the  lowest  indicated  forecast  potential 
was  assigned  a  weight  factor  of  1.0.   The  weight  factors 
for  the  other  predictors  for  each  period  were  computed  by 
dividing  the  values  of  the  associated  indicated  forecast 
potentials  by  the  lowest  potential  for  the  respective  fore- 
cast period.   This  procedure  normalized  the  numerical  weight 
factors  based  on  the  smallest  value  of  an  indicated  forecast 
potential  for  each  period.   The  selected  predictors  and 
their  respective  weight  factors  are  listed  in  Appendix  G. 

At  this  point,  the  primary  objective  of  the  study  was 
accomplished.   An  examination  of  the  potential  value  of  the 
analog  elements  as  predictors  of  intensification  was  made; 
a  number  of  predictors  was  selected  for  each  forecast  period 
and  associated  weight  factors  derived.   To  test  the  findings 
of  this  part  of  study  and  to  evaluate  the  merit  of  the 
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procedures  utilized  for  identifying  useful  predictors,  a 
forecast  scheme  was  developed. 

H.   DEVELOPMENT  OF  A  FORECAST  SCHEME 

The  forecast  scheme  to  be  developed  could  function  in 
many  different  ways.   A  weighted  average  intensification 
of  all  the  selected  analogs  could  be  computed  as  a  forecast, 
or  a  filtering  technique  could  be  applied  to  all  the  analogs 
to  determine  those  which  most  closely  resembled  the  test 
cyclone  and  apply  the  intensifications  of  these  filtered 
analogs  as  a  forecast.   A  scheme  of  the  latter  type  was  de- 
vised, as  the  actual  forecast  potentials  of  the  predictors 
would  be  more  directly  utilized.   The  following  forecast 
scheme  was  developed  and  applied  to  the  30  dependent  cases 
used  to  select  the  predictors. 

The  test  cyclone  definition  and  analog  selection  opera- 
tion was  identical  to  that  developed  for  analyzing  the  de- 
rived elements  with  one  addition.   For  each  forecast  period 
(0-24,  24-48,  48-72  hours),  the  maximum  positive  and  nega- 
tive intensifications  exhibited  by  the  selected  analog  can- 
didates were  determined. 

The  second  operation  constituted  the  actual  forecast 
scheme.   The  24-hour  forecast  was  first  determined.   All 
selected  analogs  were  classified  according  to  their  24-hour 
intensification.   The  analogs  that  intensified  between  -5 
and  +5  knots,  inclusive,  were  classed  together  and  assigned 
a  characteristic  intensification  of  zero.   Next,  all  analogs 
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that  intensified  between  +5  and  +15  knots  were  classed  to- 
gether and  assigned  a  characteristic  intensification  of 
+10  knots.   Analogs  that  exhibited  an  intensification  of 
+5  knots  were  assigned  to  both  classes  as  the  data  were 
not  accurate  enough  to  warrant  a  rigid  delineation  between 
classes.   The  procedure  was  repeated  until  all  the  inten- 
sifying analogs  were  classified.   Next,  all  the  analogs 
exhibiting  negative  intensification  were  classified  in  a 
similar  manner.   One  of  the  resulting  characteristic  inten- 
sifications would  be  selected  as  a  forecast. 

The  characteristic  intensifications  and  the  number  of 
analogs  classed  with  each  were  printed  out.   The  distribu- 
tion of  the  analogs  among  the  classes  indicated  the  most 
to  least  frequently  occurring  intensifications  of  all  past 
six-hourly  reports  of  all  past  tropical  cyclones  exhibiting 
a  time  of  occurrence,  position,  and  initial  intensity  sim-'- 
lar  to  those  of  the  test  cyclone.   The  distribution  there 
fore  represented  an  intensification  climatology  of  a 
precisely  defined  group  of  tropical  cyclone  reports.   Tht; 
characteristic  intensification  of  the  class  containing  the 
largest  number  of  analogs  was  defined  as  the  modal  intensi- 
fication . 

For  each  predictor,  the  degree  of  matching  was  deter- 
mined for  each  analog  class.   As  before,  the  degree  of  match- 
ing was  defined  as  the  number  of  actual  matches  relative  to 
the  number  of  possible  matches.   The  matching  criteria 
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(acceptance  intervals)  were  those  utilized  in  the  selection 
of  the  predictors.   The  analog  class  that  exhibited  the 
highest  degree  of  matching  for  each  predictor  was  determined, 
and  the  weight  factor  associated  with  that  predictor  was 
assigned  to  that  class  of  analogs.   If  two  or  more  classes 
exhibited  the  same  degree  of  matching  for  a  predictor,  the 
weight  factor  was  assigned  to  all  of  those  classes.   The 
assigned  weight  factors  for  each  class  were  then  summed  and 
the  summation  for  each  class  was  printed  out.   These  same 
procedures  were  then  applied  to  the  48-  and  72-hour  forecast 
period  s  . 

Originally,  it  was  the  intention  to  select  the  charac- 
teristic intensification  associated  with  the  largest  summa- 
tion of  weight  factors  as  the  forecast.   However,  analysis 
of  early  results  revealed  a  problem  area  and  a  previously 
overlooked  factor:  the  variable  probability  of  chance  ele- 
ment matching.   The  numbers  of  analogs  associated  with  the 
different  classes  varied  from  as  few  as  one  to  ten  to  as 
many  as  80  to  100,  as  an  example.   In  most  cases,  the 
classes  containing  relatively  very  few  analogs  resulted  in 
relatively  very  large  weight  factor  summations.   The  trend 
was  so  persistent  that  some  unaccounted-for  factor  was  sus- 
pected to  be  influencing  the  results.   It  was  believed 
that  the  smaller  the  number  of  analogs  in  a  given  class, 
the  greater  was  the  probability  of  getting  element  matching 
due  to  pure  chance.   The  implications  of  this  factor  are 
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discussed  further  in  the  results  section  and  only  the  pro- 
cedures attempted  to  account  for  chance  matching  are  stated 
here  . 

To  compensate  for  the  chance  matching  factor,  it  was 
necessary  to  account  for  the  number  of  analogs  in  each  class 
Simply  multiplying  the  weight  factor  summation  of  each  class 
by  the  number  of  associated  analogs  would  not  suffice,  be- 
cause the  ranges  and  orders  of  magnitude  of  the  numbers  of 
analogs  were  much  larger  than  those  of  the  weight  factor 
summations.   The  results  would  be  almost  totally  functions 
of  modal  intensification.   To  reduce  the  effects  of  the  num- 
ber of  analogs  in  each  class  on  the  results,  the  square 
roots  of  these  numbers  were  computed  and  multiplied  by  the 
respective  weight  factor  summations.   This  procedure  alle- 
viated the  problem  of  chance  matching  to  a  large  degree, 
but  in  many  cases  in  which  there  were  extremely  large  varia- 
tions in  the  numbers  of  analogs  among  the  classes,  the  prob- 
lem still  persisted.   Therefore  it  was  decided  to  eliminate, 
as  possible  forecasts,  the  characteristic  intensifications 
of  classes  with  relatively  very  few  analogs.   In  order  to 
maintain  a  strictly  objective  forecast  scheme,  a  definite 
procedure  was  needed  to  select  those  classes  whose  charac- 
teristic intensifications  would  be  retained  for  the  forecast 
determination.   Only  those  characteristic  intensifications 
whose  frequencies  of  occurrence  were  in  the  upper  50  percent 
of  the  entire  distribution  were  retained.   That  is,  the 
classes  were  selected  in  order  of  decreasing  number  of 
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associated  analogs  until  one  half  of  the  total  number  was 
selected,  unless  an  equal  number  of  analogs  in  two  or  more 
classes  required  more  than  one  half  of  the  total.   The 
development  of  this  procedure  was  purely  subjective,  but 
it  did  provide  for  an  objective  determination  of  a  forecast, 
and  a  complete  revision  of  the  forecast  scheme  in  order  to 
account  for  chance  matching  was  impractical  at  this  time. 

For  each  test  case,  three  forecasts  were  derived  based 
on  different  criteria: 

1)  modal  intensification:  characteristic 
intensification  associated  with  the 
largest  number  of  analogs 

2)  weighted  predictors:  characteristic  in- 
tensification associated  with  the  largest 
weight-factor  summation  (after  elimina- 
tion of  specified  analog  classes) 

3)  weighted  climatology:   characteristic 
intensification  associated  with  the 
largest  product  of  weight-factor  summa- 
tion and  square  root  of  the  number  of 
associated  analogs  (after  elimination 
of  specified  analog  classes). 

For  each  method  the  24-hour  intensification  was  de'.erminec 
and  added  to  the  initial  intensity  to  obtain  the  24-hour 
forecast  intensity.   Then  the  24-  to  48-hour  intensification 
was  determined  and  added  to  the  24-hour  intensity  to  obtain 
the  48-hour  forecast  intensity.   The  72-hour  intensity  was 
determined  similarly.   A  minimum  forecast  intensity  of  30 
knots  was  imposed  on  each  forecast  method.   This  was  done  to 
reduce  the  degrading  effects  of  a  highly  underforecasted  24- 
hour  intensity,  for  example,  on  the  48-  and  72-hour  forecast 
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intensities  which  might  otherwise  be  very  accurate.   The 
resulting  forecasts  were  compared  to  the  official  JTWC 
forecasts  issued  for  the  30  test  cases.   The  forecast  scheme 
was  then  applied  to  62  more  independent  test  cases  selected 
using  the  same  criteria  for  the  original  30  test  cases. 
These  test  cases  are  listed  in  Appendix  H.   The  results  of 
these  forecasts  were  also  compared  to  official  JTWC  fore- 
casts.  A  discussion  of  the  forecast  results  and  the  results 
of  the  study  in  general  is  given  in  the  next  section. 
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Ill  .   RESULTS  AND  CONCLUSIONS 

A.   DISCUSSION  OF  SELECTED  PREDICTORS 

No  attempt  was  made  to  justify  rigorously,  using  either 
meteorological  principles  or  statistical  testing,  the  analog 
elements  that  exhibited  high  indicated  forecast  potentials 
and  thus  were  selected  as  predictors.   There  are,  however, 
some  obvious  results  and  trends  worth  noting.   See  Appendix 
G  for  the  selected  predictors  for  each  forecast  period. 

Extrapolation  of  past  intensity  changes  in  operational 
forecasts  has  been  used  extensively  and  has  been  shown  to 
produce  reliable  results  for  short-term  forecasts.   Of  the 
derived  elements,  the  one  exhibiting  the  highest  indicated 
forecast  potential  for  the  24-hour  forecast  was  the  past  24- 
hour  change  in  intensity.   This  element  was  also  selected 
as  a  predictor  for  the  48-hour  forecast  period,  but  it  ex- 
hibited a  lower  indicated  forecast  potential  than  it  did 
for  the  24-hour  forecast  period.   For  the  72-hour  forecast 
period,  no  past  intensity  change  elements  exhibited  large 
enough  forecast  potentials  to  warrant  selection  as  predic- 
tors.  This  trend  agrees  with  the  fact  that  extrapolation  of 
past  intensity  changes  in  operational  forecasting  has  given 
more  accurate  results  for  the  short-term  forecasts. 

The  results  indicated  that  the  height  and  latitude  of 
the  700-mb  subtropical  ridge  at  the  longitude  of  the  cyclone 
were  related  to  intensification  during  all  forecast  periods. 
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Of  all  the  original  and  derived  elements,  these  quantities 
and  derivatives  involving  them  appeared  as  predictors  the 
most  persistently.   The  resulting  weight  factors  of  predic- 
tors involving  these  quantities  generally  increased  in  the 
most  extensive  forecast  period.   For  example,  the  weight 
factors  derived  for  the  700-mb  ridge  latitude  increased  from 
1.12  and  1.13  for  the  24-  and  48-hour  forecast  periods  to 
1.52  for  the  72-hour  forecast  period.   The  derived  analog 
element  involving  the  difference  in  700-mb  heights  of  the 
ridge  and  cyclone  divided  by  the  difference  in  latitude 
showed  a  similar  trend.   It  was  not  designated  as  a  24-hour 
predictor,  but  the  weight  factor  associated  with  this  ele- 
ment increased  from  1.13  for  the  48-hour  period  to  1.33  for 
the  72-hour  period .   The  past  24-hour  change  in  this  pre- 
dictor also  exhibited  weight  factors  increasing  with  the 
forecast  period,  e.g.,  from  1.18  to  1.37  to  1.81  in  the  ?'►-, 
48-,  and  72-hour  forecast  periods,  respectively. 

The  past  12-hour  speed  of  movement  element  displayed 
decreasing  indicated  forecast  potentials  from  the  24-  to 
72-hour  forecast  periods  and  was  defined  as  a  predictor  for 
only  the  24-  and  48-hour  periods.   The  weight  factors 
associated  with  this  predictor  decreased  from  1.42  for  the 
24-hour  forecast  to  1.07  for  the  48-hour  forecast.   The 
past  24-hour  speed  of  movement  element  was  selected  as  a 
predictor  for  the  72-hour  period,  but  possessed  the  lowest 
weight  factor  for  that  period. 
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The  700-mb  trough  height  and  the  past  24-hour  change  in 
this  quantity  dominated  the  48-hour  forecast  predictors,  in 
that  they  possessed  the  two  highest  indicated  forecast  po- 
tentials and  associated  weight  factors.  The  700-mb  trough 
height  also  possessed  the  highest  indicated  forecast  poten- 
tial of  any  of  the  original  analog  elements  for  the  72-hour 
forecast  period. 

The  estimated  1000/700-mb  thickness  element  was  not 
itself  defined  as  a  predictor  for  any  of  the  forecast  per- 
iods.  However  the  past  24-  and  48-hour  changes  of  this 
element  did  exhibit  sufficient  indicated  forecast  potentials 
to  be  designated  as  predictors.   There  was  little  observable 
trend  in  the  indicated  forecast  potential  of  these  predic- 
tors from  one  forecast  period  to  the  next.   The  past  48-hour 
change  in  thickness  was  defined  as  a  predictor  for  the  24- 
hour  period;  no  thickness  elements  were  defined  as  predic- 
tors for  the  48-hour  period;  the  past  24-hour  and  again  the 
past  48-hour  changes  in  thickness  were  designated  predictors 
for  the  72-hour  period. 

It  is  also  interesting  to  note  that  the  minimum  700-mb 
height  above  the  cyclone  was  defined  as  a  predictor  for  the 
24-hour  forecast  period  only.   However,  derived  elements 
involving  this  quantity  were  designated  as  predictors  for 
all  forecast  periods. 

The  elements  designated  as  predictors  from  the  results 
of  this  study  are  not  necessarily  considered  to  be  the 
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actual  best  predictors  of  tropical  cyclone  intensification 
among  the  possibilities  derivable  from  the  data.   An  almost 
infinite  number  of  combinations  of  the  original  18  analog 
elements  could  be  evaluated  as  possible  predictors.   As  will 
be  discussed  later,  the  results  of  just  those  elements  con- 
sidered in  this  study  could  be  altered  by  modified  and  im- 
proved procedures. 

B.   FORECAST  VERIFICATION  AND  RESULTS 

1 .   Results  of  All  Dependent  Test  Cases 

The  analog  forecasts  derived  for  the  30  dependent 
and  61  independent  test  cyclones,  as  well  as  the  available 
official  JTWC  forecasts  supplied  by  FWC  Guam  were  verified 
on  the  basis  of  two  main  criteria:  (1)  absolute  value  of  the 
forecast  intensity  errors;  (2)  the  algebraic  and  absolute 
values  of  forecast  intensity  errors  relative  to  the  verify- 
ing intensities  and  expressed  as  percent  errors.   The  first 
of  these  criteria  is  of  little  practical  value  in  an  opera- 
tional forecast  procedure.   It  considers  only  the  magnitude 
of  the  error  without  regard  to  the  actual  magnitude  of  the 
cyclone  intensity.   For  the  purposes  of  operational  forecast 
ing  and  contingency  procedures  of  all  concerned  activities, 
a  20-knot  error  in  a  24-hour  forecast,  for  example ,  would 
be  considered  misleading  and  inadequate  for  a  tropical 
storm  with  a  verifying  intensity  of  40  or  50  knots.   The 
same  error  for  the  same  forecast  period,  however,  would  be 
considered  adequate  for  a  super  typhoon  with  a  verifying 
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intensity  of  130  knots  [3].   Therefore  it  is  the  relative 
forecast  error  that  is  of  primary  concern.   The  Joint  Typhoon 
Warning  Center  has  established  relative  error  acceptability 
criteria  for  different  forecast  time  periods.   The  criteria 
for  the  24-,  48-  and  72-hour  forecast  periods  are  defined 
in  [3]  and  are  presented  in  Table  I.   The  criteria  are  re- 
laxed for  more  extended  forecast  periods  as  contingency  re- 
quirements are  more  flexible  over  longer  time  periods  [3] . 
Average  values  of  the  absolute  intensity  errors  are  provided 
only  to  indicate  to  what  extent  the  error  magnitudes  vary 
from  one  forecast  period  to  the  next  and  to  provide  an  addi- 
tional comparison  of  the  analog  forecasts  to  those  of  JTWC. 
The  analog  forecast  results  are  presented  for  all  three 
forecast  methods :modal  intensification,  intensification  based 
on  weighted  predictors,  and  intensification  based  on  weighted 
climatology  (weighted  combination  of  first  two  methods). 

Table  II  and  Table  III  present  the  average  absol  te 
forecast  errors  for  all  of  the  dependent  and  independent  test 
cases,  respectively.   Generally  for  both  set  of  cases,  t  le 
largest  increase  in  errors  from  the  24-  to  72-hour  forecast 
periods  occurred  between  the  24-  and  48-hour  forecast  per- 
iods.  Relatively  little  forecast  deterioration  was  observed 
between  the  48-  and  72-hour  periods.   This  trend  was  exhibi- 
ted by  all  forecast  methods  with  two  exceptions:  the  modal 
intensification  for  the  dependent  cases  and  the  weighted 
climatology  intensification  for  the  independent  cases. 
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Table  IV  and  Table  V  provide  the  24-hour  relative 
forecast  errors  for  the  dependent  and  independent  test 
cases.   For  the  dependent  cases  the  weighted  climatology 
produced  the  best  results,  but  for  the  independent  cases, 
modal  intensification  was  generally  the  most  accurate.  For 
the  weighted  predictors  and  weighted  climatology  intensi- 
fications, the  results  deteriorated  to  varying  degrees  in 
going  from  the  dependent  to  the  independent  test  cases. 
But  for  the  modal  intensification,  the  results  of  the  inde- 
pendent cases  showed  improvement  over  those  of  the  dependent 
cases,  with  the  exception  of  the  10%  error  category. 

For  both  the  dependent  and  independent  test  cases, 
the  weighted  predictors  intensification  showed  the  greatest 
fraction  of  useful  forecasts  (_<  30%  error)  that  were  accu- 
rate to  the  measurement  of  the  data  (<_  10%  error)  .   For 
example,  considering  the  dependent  cases,  (Table  V),  this 
forecast  method  resulted  in  7 5% ( 50 / 6 7= . 75 )  of  the  useful 
forecasts  being  accurate  to  the  measurement  of  the  data. 
These  percentages  for  the  modal  and  weighted-climatology 
intensifications  for  the  dependent  cases  were  55%  and  62%, 
respectively.   This  indicates  that  although  the  weighted- 
predictors  intensification  did  not  produce  as  many  useful 
forecasts  as  the  other  methods,  it  became  relatively  more 
skillful  than  the  other  methods  as  the  forecast  error 
decreased . 

The  48-hour  relative  forecast  errors  are  presented 
in  Table  VI  and  Table  VII.   The  modal  intensification  and 
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that  of  weighted  climatology  exhibited  nearly  the  same  re- 
sults, especially  for  the  useful  (<_  40%  error)  and  adequate 
(<_  30%  error)  forecast  categories.   All  forecast  methods 
suffered  deterioration  in  going  from  the  dependent  to  inde- 
pendent test  cases.   The  modal  intensification  exhibited 
the  largest  fraction  of  useful  forecasts  that  were  accurate 
to  the  measurement  of  the  data  for  both  sets  of  test  cases. 
The  weighted  predictors  showed  the  smallest  fraction  for 
the  dependent  cases  while  the  weighted  climatology  showed 
the  smallest  fraction  for  the  independent  cases. 

The  relative  forecast  errors  for  the  72-hour  fore- 
cast period  are  presented  in  Table  VIII  and  Table  IX.  For 
the  dependent  cases,  the  best  results  were  observed  for 
the  modal  and  weighted  climatology  intensifications,  while 
modal  intensification  dominated  with  the  independent  tests. 
Again  the  results  generally  deteriorated  from  the  dependent 
to  independent  cases  for  all  forecast  methods,  but  to  a 
lesser  extent  for  the  modal  intensification.   The  modal 
intensification  again  exhibited  the  largest  fraction  of 
useful  forecasts  that  were  accurate  to  the  measurement  of 
the  data  for  both  sets  of  test  cases.   As  with  the  48-hour 
forecasts,  the  weighted  predictors  showed  the  smallest 
fraction  for  the  dependent  cases,  and  the  weighted  clima- 
tology showed  the  smallest  fraction  for  the  independent 
cases.   The  fraction  associated  with  the  weighted  predictors 
showdd  an  increase  from  the  dependent  to  independent  test 
cases  (16/64  =  .25  to  16/55  =  .29). 
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Table  X  presents  the  relative  forecast  errors  of  the 
independent  test  cases  stratified  by  initial  intensity.   Two 
error  categories  are  presented  —  those  that  resulted  in  in- 
adequate forecasts,  and  those  that  resulted  in  forecasts 
accurate  to  the  measurement  of  the  data.   The  test  cases 
were  stratified  on  the  basis  of  whether  or  not  their  initial 
intensities  were  of  typhoon  strength.   For  the  24-  and  48- 
hour  forecast  periods  the  number  of  occurrences  of  inadequate 
forecasts  was  substantially  less  for  the  typhoon  cases  than 
it  was  for  the  storm  or  depression  cases  for  each  of  the 
three  forecast  methods.   But  for  the  72-hour  forecast  period 
the  results  were  much  more  similar  for  the  modal  intensifi- 
cation, and  essentially  identical  for  the  remaining  two 
methods.   The  number  of  forecasts  accurate  to  within  data 
measurements  for  each  forecast  period  was  generally  greater 
for  the  typhoon  cases  than  for  the  storm  or  depression 
cases.   For  the  24-hour  forecast  period  the  number  of  for  - 
casts  accurate  to  the  measurement  of  the  data  was  signifi 
cantly  greater  for  the  typhoon  cases  than  for  the  depression 
or  storm  cases  for  each  forecast  method.   Except  for  the 
modal  intensification,  this  stage-related  trend  decreased 
substantially  in  the  more  extended  forecast  periods.   These 
results  imply  that  the  intensifications  of  the  more  intense 
tropical  cyclones  are  generally  more  predictable  than  those 
of  the  weaker,  developing  or  dying,  depressions  and  storms. 
The  results  in  Table  X  are  based  on  relative  errors,  and 
therefore  do  not  imply  that  the  absolute  forecast  errors 
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decrease  with  increasing  initial  intensity.   But  as  stated 
previously,  an  absolute  error  of  20  knots  is  not  nearly  as 
significant  for  a  130-knot  typhoon  as  it  is  for  a  40-  to 
50--knot  storm.   Therefore,  if  it  is  accepted  that  only  rela- 
tive errors  are  meaningful,  the  implications  of  the  results 
in  Table  X  are  valid. 

Frequency  histograms  of  relative  errors  for  the 
three  forecast  methods  were  constructed  for  both  dependent 
and  independent  test  cases.   The  histograms  for  the  weighted' 
climatology  forecast  method  are  representative  of  all  those 
constructed  and  are  presented  in  Figures  1  through  6.   In 
all  cases,  relative  errors  of  +  5%  or  less  were  classified 
as  zero  forecast  error.   For  positive  errors,  N%  error  in- 
cluded all  errors  within  (N-4)%  and  (N+5)%  inclusive,  and 
for  negative  errors,  -N%  error  included  all  errors  within 
(-N+4)%  and  (-N-5)%,  inclusive.   Therefore  the  results 
presented  in  the  histograms  do  not  exactly  agree  with  those 
presented  in  the  tables. 

The  histograms  for  the  dependent  cases  reveal  that 
the  forecast  ability  of  the  weighted-climatology  method 
decreased  from  the  24-  to  48-  to  72-hour  forecast  periods. 
That  is,  the  numbers  of  occurrences  of  small  relative  fore- 
cast errors  (<_  20%)  decreased  and  the  distribution  became 
more  variable  for  the  more  extended  forecast  periods.   No 
persistent  skewness  toward  positive  or  negative  errors 
resulted.   However  it  was  observed  that  the  largest  forecast 
errors  were  consistently  positive  errors. 
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The  results  for  the  independent  test  cases  differed 
somewhat  from  the  dependent  cases.   The  most  obvious  differ- 
ence was  the  persistent  skewness  toward  negative  forecast 
errors,  especially  in  the  24-  and  48-hour  forecast  periods. 
The  fact  that  this  result  did  occur  might  have  made  possible 
the  use  of  an  adjustment  factor  or  function  to  modify  and 
improve  the  results.   This  was  not  attempted,  however,  in 
this  study.   The  histograms  reveal  a  decrease  in  the  number 
of  occurrences  of  small  forecast  error  with  increasing  fore- 
cast periods.   One  similarity  of  these  results  to  those  of 
the  dependent  cases  was  the  fact  that  the  largest  errors 
were  consistently  positive  errors. 

2  .   Comparison  of  Analog  Forecasts  to  Official 
JTWC  Forecasts 

As  a  further  measure  of  the  significance  of  the 
forecast  results,  comparisons  of  the  analog  forecasts  to 
those  of  the  official  JTWC  forecasts  were  made  for  homo- 
geneous sets  of  dependent  and  independent  test  cases.   Com- 
parisons were  made  for  the  24-  and  48-hour  forecast  periods 
only  due  to  the  very  limited  number  of  applicable  72-hour 
forecasts  available  from  FWC  Guam.   The  analog  results  were 
recomputed  based  on  only  those  test  cases  for  which  the 
official  forecasts  were  available. 

The  average  absolute  forecast  errors  are  presented 
in  Table  XI  and  Table  XII  for  the  dependent  and  independent 
test s , respect ively .   As  the  figures  indicate,  the  average 
errors  resulting  from  all  three  analog  methods  were  from  35% 


53 


less  than,  to  at  worst  equal  to  those  of  the  official  fore- 
casts for  the  dependent  test  cases.   For  the  independent 
cases,  the  errors  for  the  analog  forecasts  were  from  17%  to 
almost  60%  greater  than  official  forecast  errors.   The 
weighted  predictors  method  gave  the  worst  results  for  the 
independent  cases.   Neglecting  this  method,  the  largest  per- 
cent increase  in  error  was  33%  (for  the  weighted  climato- 
logy) .   The  overall  best  results,  considering  both  forecast 
periods  and  both  sets  of  test  cases,  were  produced  by  the 
modal  intensification.   However,  the  weighted  predictors 
and  weighted  climatology  did  result  in  smaller  average  fore- 
cast errors  for  the  24-hour  dependent  forecasts  than  did 
the  modal  intensification. 

Table  XIII  shows  the  comparisons  of  relative  fore- 
cast errors  among  the  analog  methods  and  the  JTWC  forecasts 
for  the  dependent  tests.   These  results  reveal  that  all  of 
the  analog  methods  exhibited  some  equal  or  superior  foreca  t 
ability  over  the  official  forecasts.   Exceptions  were  the 
modal  and  weighted-predictors  24-hour  forecasts  for  the 
<_  30%  error  category,  and  the  weighted-predictors  48-hour 
forecasts  for  the  <_  20%  and  <_  10%  error  categories.   All  but 
one  of  the  analog  methods  exhibited  a  greatest  fraction  of 
useful  forecasts  that  were  accurate  to  the  measurement  of 
the  data  for  both  periods  than  did  the  official  forecast. 
This  indicates  that  as  the  forecast  error  decreased,  the 
analog  methods  (with  one  exception)  became  relatively  more 
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skillful  than  did  the  official  forecasts.   The  one  exception 
to  this  result  was  the  weighted  predictors  intensification 
for  the  48-hour  period  only.   However,  the  fraction  exhibited 
by  the  weighted  predictors  method  for  the  24-hour  forecast 
period  was  substantially  greater  than  any  of  the  others. 

The  relative  error  comparisons  for  the  independent 
test  are  shown  in  Table  XIV.   For  the  most  part,  the  offi- 
cial forecasts  proved  to  be  superior  to  the  analog  fore- 
casts.  Exceptions  were  the  24-hour  weighted  climatology 
forecasts  for  the  _<  10%  category  error  category  and  the 
48-hour  modal  and  weighted  climatology  forecasts  for  the 
'<    10%  error  category.   For  all  analog  methods  for  both  fore- 
cast periods,  the  superiority  of  the  official  forecasts 
decreased  with  decreasing  error  categories  from  <_    30%  to 
j<  10%  errors.   However,  the  trend  is  just  the  opposite  in 
going  from  the  _<  40%  to  <^  30%  error  categories.   With  one 
exception,  the  fractions  of  the  useful  forecasts  that  were 
accurate  to  the  measurement  of  the  data  were  greater  for 
the  analog  forecasts  than  for  the  official  forecasts. 
Again,  the  exception  was  the  weighted  predictor  method  for 
the  48-hour  forecast  period. 

C.   CONCLUSIONS 

1 .   Implications  of  Forecast  Results 

The  overall  results  of  the  forecast  verification 
are  considered  competitive  with  official  forecasts  and  thus 
encouraging.   The  general  deterioration  of  the  results  from 
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the  dependent  to  independent  test  cases  was  expected;  none- 
theless substantial  numbers  of  useful  and  adequate  analog 
forecasts  were  generated  for  these  independent  cases.  The 
results  were  especially  encouraging  for  the  48-  and  72-hour 
forecasts,  as  these  extended  periods  offer  the  forecaster 
the  greatest  challenge.   If  these  results  were  maintained 
under  conditions  of  operational  forecasting,  they  would  at 
least  provide  the  forecaster  with  a  starting  point,  other 
than  linear  extrapolation  and  cl imat ology , f or  deriving  his 
prediction.   With  a  somewhat  reliable  starting  point,  the 
finished  forecast  might  well  be  more  accurate. 

Of  the  three  analog  forecast  methods  applied,  those 
involving  the  use  of  various  analog  elements  possess  the 
greatest  potentials  for  improvement.   Modal  intensification 
might  be  improved  slightly  by  experimenting  with  the  screen- 
ing criteria,  but  little  else  could  be  done  to  reduce  the 
forecast  errors.   However,  an  almost  unlimited  amount  of 
experimentation  with  the  analog  elements  and  their  forecast 
potentials  is  possible. 

It  has  been  established  that  the  analog  forecast 
methods  performed  more  accurately  on  test  cases  used  to 
derive  the  predictors  and  their  associated  weight  factors. 
The  analog  elements  associated  with  the  small  sample  of 
test  cases  utilized  for  this  purpose  undoubtedly  are  not 
representative  of  all  the  element  distributions.   However, 
if  all  available  past  data  were  utilized  in  deriving  the 
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predictors  and  their  weight  factors,  the  resulting  distri- 
butions of  the  analog  elements  would  be  very  representative 
of  the  possible  values  that  each  element  could  assume. 
Therefore,  the  likelihood  would  be  high  that  all  but  extreme- 
ly uncommon  cyclones  (for  which  forecasts  must  be  issued), 
would  resemble  at  least  some  cyclones  used  to  derive  the 
predictors  and  weight  factors.   Consequently,  the  relation- 
ship of  any  given  current  cyclone  to  the  information  used 
in  deriving  the  forecast  scheme  would  approach  the  relation- 
ship between  the  dependent  test  cases  and  the  information 
used  in  deriving  their  forecasts.   Thus  it  is  strongly 
believed  that  more  accurate  results  could  be  obtained  by 
increasing  to  the  maximum  extent  the  number  of  test  cases 
utilized  in  deriving  the  predictors  and  their  weight  factors 

As  stated  earlier,  the  number  of  elements  derivable 
from  the  original  18  analog  elements  is  enormous.   Evalua- 
tion of  elements  other  than  those  considered  in  this  stuc  ' 
might  well  lead  to  potentially  better  predictors  and  for*  - 
cast  results.   Also  the  technical  procedures  of  element 
evaluation  and  forecasting  could  be  modified,  or  completely 
different  avenues  of  approaches  followed. 

2  .   Evaluation  of  Investigative  Procedures 

The  encouraging  results  of  the  analog  forecasts 
derived  in  this  study  strongly  suggest  that  the  methods 
utilized  in  defining  the  predictors  and  computing  associated 
weight  factors  possess  at  least  some  degree,  of  skill.   In 
spite  of  some  crude  and  intuitively  derived  procedures 
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involved  in  the  method,  the  resulting  predictors  apparently 
do  possess  at  least  some  degree  of  forecast  potential. 
Examples  of  these  crude  procedures  are  the  definitions  of 
the  acceptance  intervals  for  original  and  derived  element 
matching.   A  "sledge  hammer"  effort  was  also  applied  in  the 
forecast  method  to  alleviate  the  problems  of  variable 
chance  element  matching. 

The  most  severe  shortcoming  of  the  methods  utilized 
was  the  failure  to  properly  account  for  the  variable  proba- 
bility of  chance  element  matching.   The  effects  of  this 
factor  were  not  realized  until  forecast  outputs  were 
analyzed.   However,  it  is  likely  that  the  same  distorting 
effects  were  present  in  the  element  investigation  procedures 

In  the  second  and  third  approaches  to  the  main  ob- 
jective, the  selected  analogs  were  classified  on  the  basis 
of  a  maximum  allowable  deviation  in  intensification  from 
that  of  the  test  case.   First,  all  analogs  that  intensified 
within  +  10  knots  per  24-hour  period  were  classified  and 
the  degree  of  element  matching  determined.   Then  all  analogs 
that  underwent  a  change  in  intensity  of  +  20  knots,  +  30 
knots  and  finally  +  40  knots  per  24-hour  period  were  classi- 
fied and  for  each  maximum  allowable  intensification  devia- 
tion, the  degrees  of  element  matching  were  computed.   It  is 
obvious  that  as  these  deviations  increased,  the  number  of 
acceptable  analogs  increased,  which  in  turn,  generally  in- 
creased the  number  of  possible  analog  element  matches. 
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For  a  particular  analog  element,  a  population  could 
be  defined  as  the  entire  set  of  values  for  that  element  coil' 
tained  in  the  complete  data  file.   The  different  classes  of 
analogs  represent  samples  of  various  sizes  of  this  element 
population.   Generally  as  the  number  of  analogs  increased, 
with  increasing  intensification  range,  the  size  of  the  ele- 
ment samples  also  increased.   The  variable  sample  sizes  is 
one  factor  of  the  chance  element  matching. 

The  other  factor  is  the  location,  within  an  element 
population,  of  the  value  of  the  test  case  analog  element. 
That  is,  in  general,  the  farther  (closer)  the  test  case 
value  is  from  the  expected  (i.e.  mean)  value,  the  smaller 
(greater)  the  probability  of  getting  a_  match  by  pure  chance 

In  the  procedures  of  this  study,  the  degree  of 
matching  was  defined  as  the  percent  of  matches  to  the  total 
number  of  possible  matches.  However,  the  ratio  involved  in 
this  computation  can  also  be  considered  an  average  of  ranks 
This  ratio  is  nothing  more  than  the  sum  of  ranks  (0  or  1) 
divided  by  the  number  of  ranks. 

The  Central  Limit  Theorem  can  be  applied  to  this 
problem.   The  implications  are  that  as  the  sample  size  in- 
creases, and  as  the  test  cyclone  element  value  departs  the 
expected  value  in  the  element  population,  the  smaller  the 
probability  that  a  specified  average  rank  (e.g.  .50)  will 
occur  due  to  chance  alone.   This  implies  a  smaller  proba- 
bility of  a  specified  degree  of  matching  (e.g.  50%)  due  to 
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chance  as  the  sample  size  increases.   But  as  the  sample  size 
increases  and  as  the  test  case  element  value  approaches  the 
expected  value  in  its  distribution,  the  greater  the  proba- 
bility of  a  specified  average  rank  due  to  chance.   This 
implies  a  greater  probability  of  a  specified  degree  of  match- 
ing,due  to  chance  alone,  as  the  sample  size  increases.   The 
significance  of  a  match  varies  with  the  location  of  the  test 
case  element  value  in  its  distribution.   To  account  for  this, 
the  size  of  the  acceptance  interval  should  vary  from  smaller 
magnitudes  in  high  density  regions  of  the  population  to 
larger  magnitudes  in  low  density  regions. 

To  properly  compensate  for  chance  element  matching, 
both  the  analog  class  sizes  and  the  location  of  the  parti- 
cular test  case  element  value  within  its  population  must  be 
accounted  for.   This  requires  a  thorough  study  of  the  dis- 
tributions of  all  predictors  to  be  utilized. 

Since  chance  matching  could  tend  to  increase  the 
degree  of  matching  for  both  large  and  small  analog  classes 
(depending  on  the  location  of  the  test  case  element  value 
in  its  population),  the  procedures  taken  in  this  study  to 
compensate  for  chance  matching  were  only  partially  valid. 
The  procedures  were  proper  only  when  the  test  case  element 
values  were  in  low  density  regions  of  their  distribution. 

Regardless  of  the  shortcomings  of  the  procedures  in 
this  study,  the  results  warrant  further  study   into  the 
potential  of  analog  schemes  in  forecasting  tropical  cyclone 
intensification.   It  is  felt  that  concentrated  efforts  in 
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this  area  could  lead  to  improved  and  more  significant  re- 
sults, and  the  evolution  of  a  reliable  and  efficient  in- 
tensification forecast  aid. 

3 .   Thoughts  on  Subsequent  Studies 

In  all  the  investigative  and  forecast  procedures 
utilized  in  this  study,  the  first  operation  was  to  screen 
out  those  six-hourly  reports  of  past  tropical  cyclones  that 
exhibited  similarity  to  the  test  case  with  respect  to  loca- 
tion and  time  of  occurrence.   Further  screening  criteria 
then  determined  those  analog  candidates  that  resembled  the 
test  case  with  respect  to  the  analog  elements.   The  ideal 
analog  procedure  would  be  to  continue  the  screening  criteria 
until  one  analog  candidate  was  defined  to  be  the  most  simi- 
lar to  the  test  case.   This  ideal  procedure  would  require 
that  the  data  observations  be  exact.   It  would  also  be 
necessary  to  know  exactly  what  parameters  are  the  ideally 
best  predictors.   These  ideally  best  predictors  would  have 
to  possess  precise  relationships  to  intensification,  and 
the  relationships  would  have  to  hold  true  for  all  cyclone 
occurrences.   In  reality  the  data  are  not  exact,  and  reliable 
predictors  with  precise  relationships  to  intensification 
have  not  yet  been  defined.   This  forces  the  analog  technique 
to  rely  on  average  behavior  in  order  to  generate  a  quanti- 
tative forecast.   With  the  use  of  averages,  it  is  hoped 
that  random  observational  errors  and  non-persistent   rela- 
tionships will  be  smoothed  to  a  degree  that  enables  the 
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utilization  of  the  data  in  an  adequate  forecast  scheme. 
Indeed,  if  the  data  were  exact  and  precise  relationships 
between  predictors  and  intensification  were  known,  deter- 
ministic techniques  would  replace  statistical  schemes. 
To  generate  a  forecast,  it  would  be  necessary  to  simply 
substitute  predictor  values  in  an  equation  and  solve  for 
the  intensification. 

It  is  felt,  however,  that  some  efforts  should  be 
made  to  strive  for  a  more  ideal  analog  forecast  technique, 
involving  more  sophisticated  screening  criteria  aimed  at 
reducing  the  number  of  analog  cyclones  used  in  generating 
the  forecast.   A  continuous- type  degree  of  matching 
criterion  was  attempted  in  the  first  approach  of  this  study 
Due  to  the  inconclusive  results  of  this  procedure,  a  match, 
no-match  criterion  was  investigated.   A  combination  of 
these  two  matching  criteria  might  lead  to  better  results. 
That  is  for  each  analog  element  apply  a  match,  no-match 
screening  procedure  to  eliminate  grossly  unmatched  values, 
and  to  identify  those  values  within  some  prescribed  range 
from  the  test  case  value.   With  a  large  degree  of  noise 
eliminated,  then  proceed  with  a  cont inuous -type  degree  of 
matching  operation.   For  example,  the  difference  in  element 
values  relative  to  the  magnitude  of  the  acceptance  interval 
(used  in  the  first  step)  could  be  examined.   This  type  of 
procedure  could  possibly  define  more  precise  and  persistent 
relationships  between  analog  elements  and  intensification 
than  those  which  were  defined  in  this  study.   With  the  use 
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of  more  precise  relationships,  the  screening  procedure  in 
the  forecast  scheme  could  be  more  restrictive  and  thereby 
select  a  fewer  number  of  analogs  that  are  more  similar  to 
the  test  case  with  respect  to  predictor  values.   In  any 
event,  further  studies  might  reveal  to  what  extent  average 
behavior  must  be  relied  upon  to  smooth  out  the  data  errors 
and  non-precise  and  quasi-persistent  relationships. 
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TABLE  I 


Forecast-intensity  relative  error  acceptability  criteria, 
as  defined  by  JTWC. 


E  =  forecast-intensity  error  relative 
to  verification  intensity 


Forecast  Value 


Forecast  Period 


24-hour 


48-  and  72-hour 


Accurate  to  within 
intensity  measure- 
ment error 


E  <  10% 


E  <  10% 


Ad  equat  e 


E  <  20% 


E  <  30% 


Useful 


E  <  30% 


E  <  40% 


Inadequat  e 


E  >  3  0% 


E  >  40% 
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TABLE  II 


Average  absolute  forecast-intensity  errors  (kt)  for  all  de- 
pendent test  cyclones.  Numbers  in  parentheses  indicate  the 
number  of  applicable  forecasts  for  each  forecast  period. 


Forecas  t 

Period 

(Hr) 


Modal 


Forecast  Method 


Weigh  ted 
Predictor  s 


Weighted 
Climatology 


24 

(30) 

16 

48 

(29) 

16 

72 

(25) 

19 

13 
23 
23 


13 
20 
21 


TABLE  III 


Average  absolute  forecast-intensity  errors  (kt)  for  all  in- 
dependent test  cyclones.   Numbers  in  parentheses  indicate 
the  number  of  applicable  forecasts  for  each  forecast  period 


Forecast 

Period 

(Hr) 


Modal 


Forecast  Method 

Weighted      Weighted 
Predictors    Climatology 


24  (61)  16 
48  (59)  26 
72   (51)       28 


19 
28 
33 


26 
28 
32 
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TABLE  IV 


Percent  of  all  dependent  test  cyclones  exhibiting  the  indi' 
cated  relative  error  (E)  in  the  24-hour  forecast  intensity 
Results  based  on  30  forecasts. 

E  =  forecast-intensity  error  relative  to  verification 
intensity . 

Forecast 
Method  E(%) 


<  40  £  30  £  20  <_   10 

Modal  87  73  57  40 

Weighted 

Predictors  83  67  63  50 

Weight  ed 

Climatology  90  80  70  50 


TABLE  V 


Percent  of  all  independent  test  cyclones  exhibiting  the 
indicated  relative  errors  (E)  in  the  24-hour  forecast 
intensity.   Results  based  on  61  forecasts. 

E  =  forecast-intensity  error  relative  to  verification 
intensity . 


Forecast 
Method  E(%) 

<  40    <  30    <  20    <  10 


Modal  90      82      64      31 

Weighted 

Predictors  79      66      54      38 

Weight  ed 

Climatology  84      74      54      36 
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TABLE  VI 

Percent  of  all  dependent  tests  cyclones  exhibiting  the  indi- 
cated relative  errors  (E)  in  the  48-hour  forecast  intensity 
Results  based  on  29  forecasts. 

E  =  f orecas t -int ens ity  error  relative  to  verification 
intensity . 


Forecast 
Method  E(%) 

<  40    £  30    <  20    <_  10 

Modal  79      70      59      41 

Weigh  ted 

Predictors  70      59      43      21 

Weighted 

Climatology  76      70      59      34 


TABLE  VII 

Percent  of  all  independent  test  cyclones  exhibiting  the 
indicated  relative  errors  (E)  in  the  48-hour  forecast 
intensity.   Results  based  on  59  forecasts. 
E  ■  forecast-intensity  error  relative  to  verification 
int ens  it y . 


Forecast 
Method  E(%) 

1   A0    1  30    1  20    £  10 

Modal  71      54      44      23 

Weighted 

Predictors  64      49      30      20 

Weighted 

Climatology  68      54      34      15 
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TABLE  VIII 

Percent  of  dependent  test  cyclones  exhibiting  the  indicated 
relative  errors  (E)  in  the  72-hour  forecast  intensity. 
Results  based  on  25  forecasts. 

E  =  forecast-intensity  error  relative  to  verification 
intensity . 

Forecast 
Method  E(%) 


<  40    <  30    <  20    <  10 


Modal  68      56      52      40 

Weigh  ted 

Predictors  64      56      32      16 

Weigh  ted 

Climatology  80      56      36      24 


TABLE  IX 

Percent  of  all  independent  test  cyclones  exhibiting  the 
indicated  relative  errors  (E)  in  the  72-hour  forecast 
intensity.   Results  based  on  51  forecasts. 

E  =  forecast-intensity  error  relative  to  verification 
intensity . 


Forecast 
Method  E(%) 

<  40    <  30    <  20    <  10 


Modal  67      53      43      24 

Weighted 

Predictors  55      41      33      16 

Weighted 

Climatology  59      35      31      16 
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Figure  1.   Frequency  historgram  of  24-hour  relative  fore- 
cast errors  for  weighted  climatology  forecast  method, 
dependent  test  cyclones. 
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Figure  2.   Frequency  histogram  of  48-hour  relative  forecast 
errors  for  weighted  climatology  forecast  method,  dependent 
test  cyclones . 
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Figure  3.   Frequency  histogram  of  72-hour  relative  forecast 
errors  for  weighted  climatology  forecast  method,  dependent 
test  cyclones . 
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Figure  4.   Frequency  histogram  of  24-hour  relative  forecast 
errors  for  weighted  climatology  forecast  method,  independent 
test  cyclones . 
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Figure  5.   Frequency  histogram  of  48-hour  relative  forecast 
errors  for  weighted  climatology  forecast  method,  independent 
test  cyclones . 
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Figure  6.  Frequency  histogram  of  72-hour  relative  forecast 
errors  for  weighted  climatology  forecast  method,  independent 
test  cyclones . 
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APPENDIX  A 

Original  Analog  Parameters  and  Error 
Estimates  of  the  Analog  Elements 


The  following  parameters  were  compiled  at  six-hourly 
intervals  for  use  with  the  TYFOON  analog  forecast  scheme. 
Estimated  maximum  observational  errors  in  the  analog 
elements  were  provided  through  informal  discussion  with 
Commander  J.  D.  Jarrell,  USN.   Commander  Jarrell  is  a  co- 
originator  of  the  TYFOON  analog  forecast  scheme  and  has 
served  as  Operations  Officer  for  the  Joint  Typhoon  Warning 
Center,  Guam.   Error  estimates  for  the  past  movement  ele- 
ments are  based  on  estimated  errors  of  six  nautical  miles 
in  the  latitude  and  longitude  of  a  cyclone  position. 
Identification  Elements 

1.  Storm  name  (a  number  is  given  unnamed  storms) 

2.  Storm  identity  number.   Month,  year  and  cyclone 
sequence  number 

3.  Observation  number 
A.   Year 

5.  Month/Day 

6.  Hour  (GMT) 

7 .  Latitude 

8 .  Longitude 
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Analog  Elements 

1.  Past  12-hr  direction 
of  movement 

2.  Past  12-hr  speed 
of  movement 


Estimated  Maximum 
Observational  Errors 


6  degrees 


1  knot 


3.  Past  18-hr  direction 
of  movement 

4.  Past  18-hr  speed  of 
movemen  t 


<  6  degrees 


<  1  knot 


5.  Past  24-hr  direction 
of  movement 

6.  Past  24-hr  speed  of 
movement 


<  6  degrees 


<  1  knot 


7.  Past  48-hr  direction 
of  movement 

8.  Past  48-hr  direction  of 
movement 


<  6  degrees 


<  1  knot 


9.   Size  (average  radius  of 
outer  closed  isobar  in 
whole  degrees  of  latitude) 

10.  Past  12-hr  change  of  size 
in  degrees  latitude 

11.  Minimum  observed  sea-level 
pressure  in  whole  millibars 

12.  Past  12-hr  change  in  minimum 
sea-level  pressure 

13.  Maximum  wind 


2  degrees  latitude 


4  degrees  latitud 


5  millibars 


10  millibars 

10%  of  the  observed 
value 


14.   Minimum  700-mb  height 
(above  the  cyclone) 

15    Latitude  of  the  700-mb 

ridge  north  of  the  cyclone 

16.   700-mb  height  at  the  ridge 
line  north  of  the  cyclone 


10  meters 


2  degrees  latitude 


20  meters 


81 


17.   Longitude  at  35N  of  the 
nearest  700-mb  trough 
west  of  the  cyclone 


2  degrees  longitude 


7  00-mb  height  at  the 
intersection  of  the  trough 
line  and  35N 


30  meters 
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APPENDIX  B 
Dependent  Test  Cyclones 

Dependent  test  cyclones.   The  following  tropical  cyclone 
reports  were  selected  to  examine  the  forecast  potential  of 
the  analog  elements. 


Name 


Betty 

Cora 

Doris 

Wendy 

Wendy 

Shirley 

Shirley 

Agnes 

Agnes 

Agnes 

Agnes 

Agnes 

Agnes 

Agnes 

Polly 

Polly 

Polly 

Polly 

Louise 

Louise 

Marge 

Opal 

Opal 

Opal 

Rita 

Rita 

Carmen 

Carmen 

Carmen 

Carmen 


MO-DA-YR 


8-06-69 
8-17-69 
8-02-69 
8-29-68 
8-31-68 
8-17-68 
8-18-68 
8-30-68 
9-01-68 
9-03-68 
9-04-68 
9-05-68 
9-07-68 
9-08-68 
8-07-68 
8-09-68 
8-11-68 
8-13-68 
8-20-67 
8-21-67 
8-26-67 
9-01-67 
9-03-67 
9-11-67 
8-04-66 
8-05-66 
8-10-63 
8-12-63 
8-13-63 
8-14-63 


Time 

(GMT) 

1200 

0000 

0600 

1800 

0000 

1800 

1200 

0600 

0000 

0600 

0600 

0600 

0600 

0600 

1200 

0000 

0000 

0000 

0000 

0000 

1200 

0600 

1200 

1200 

1200 

1200 

0000 

1800 

0600 

0000 

Initial 

Int ens  i ty 

(kt) 

45 

50 

30 

100 

140 

30 

50 

45 

85 

130 

145 

115 

70 

60 

40 

55 

50 

55 

30 

55 

100 

75 

155 

90 

60 

75 

40 

125 

110 

85 
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APPENDIX  C 

Acceptance  Intervals  Defined  for  the 
Original  Analog  Elements 


The  acceptance  intervals  defined  for  the  18  original 
analog  elements  are  listed  below.   The  range  of  each  in- 
terval is  defined  as  plus  or  minus  one-fourth  of  one 
standard  deviation  of  the  respective  element's  distribu- 
tion over  the  entire  data  file  obtained  for  the  study. 
This  range  is  centered  on  the  test  cyclone  element  value. 
If  this  fraction  of  the  standard  deviation  resulted  in  a 
number  less  than  one,  the  acceptance  interval  was  defined 
as  +  1.0.   Intervals  vice  discrete  values  were  preferred 
in  order  to  prevent  the  necessity  of  identical  element 
values  in  order  for  a  match  to  occur.   Refer  to  Appendix 
A  for  analog  element  definition. 


Analog  Element  Number 

1 

2 

3 

4 

5 

6 

7 

8 

9 
10 
11 
12 
13 
14 
15 
16 
17 
18 


Acceptance  Interval 

+ 

7  .1 

degrees 

+ 

1.4 

knot  s 

+ 

7  .6 

degrees 

+ 

1.3 

kno  t  s 

+ 

7.6 

degrees 

+ 

1.2 

knot  s 

+ 

7.3 

degrees 

+ 

1.0 

knots 

+ 

1.0 

degrees 

latitude 

+ 

1.0 

degrees 

lat  itude 

+ 

5.7 

millibars 

+ 

2.1 

millibars 

+ 

7.9 

kno  t  s 

+ 

4.7 

d  ekameter s 

+ 

1.4 

degrees 

latitude 

+ 

1.0 

dekame t  < 

ars 

+ 

3.5 

degrees 

longitude 

+ 

1.0 

dekame t  < 

srs 
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APPENDIX  D 

Indicated  Forecast  Potentials  of  the 
Original  Analog  Elements 


Indicated  forecast  potentials  of  the  original  analog 
elements  for  each  forecast  period  are  listed  below.   The 
indicated  forecast  potentials  are  expressed  as  a  percent 
of  the  total  number  of  possible  occurrences  in  which  the 
trend  of  decreasing  degree  of  matching  with  increasing 
forecast  error  occurred  for  each  element.   The  results  are 
based  on  thirty  24-hour  forecasts,  twenty-nine  48-hour  fore^ 
casts  and  twenty-five  72-hour  forecasts.   Refer  to  Appendix 
A  for  element  definitions. 


Analog  Element 
Numb  er 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 


Forecast  Potential  for 
Given  'Forecast  Periods  (%) 


24-hr 

17 
41 
17 
21 
25 
21 
24 
12 
33 
13 
23 
47 
27 
33 
37 
40 
33 
20 


48-hr 

29 
32 
29 
18 
30 
19 
29 
17 
28 
28 
31 
21 
38 
31 
34 
38 
28 
41 


72-hr 

8 
12 
12 

8 
21 
21 
14 

9 
32 
28 
16 
24 
20 
20 
12 
32 
20 
36 
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APPENDIX  E 

Derived  Analog  Elements  and  their 
Acceptance  Intervals 


The  derived  analog  elements,  formed  by  various  combi- 
nations of  some  of  the  original  analog  elements,  are  listed 
below  together  with  the  acceptance  interval  tStat  was 
selected  for  each.   The  acceptance  interval  for  each  ele- 
ment is  centered  on  the  test  cyclone  value  of  that  element. 


Derived  Analog 
Element  s 

1.  (700-mb  ridge  height-minimum 

700-mb  height) 

2.  (700-mb  ridge  lat-cyclone  lat) 


Acceptance 
Intervals 


+  10  m 


+  2    deg  lat 


3. 


A. 


6. 

7. 
8. 

9. 


700-mb  ridge  height-minimum  700-mb  height 
700-mb  ridge  lat-cyclone  lat 


(700-mb  trough  height-minimum 
700-mb  height) 


5.   (700-mb  trough  long-cyclone  long) 


+  10  m/deg  lat 

+  10  m 

+  2  deg  long 


700-mb  trough  height-minimum  700-mb  height 
700-mb  trough  long-cyclone  long 


1000/700-mb  thickness 

Past  24-hr  change  in  derived 
element  #2 

Past  24-hr  change  in  derived 
element  #3 


+  10  m/deg  long 
+   2  deg  lat 

+   2  deg  lat 

+  10  m/deg  lat 
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10.  Past  24-hr  change 
element  //5 

11.  Past  24-hr  change 
element  //  6 

12.  Past  24-hr  change 
element  #7 

13.  Past  24-hr  change 
(maximum  wind) 

14.  Past  24-hr  change 
pressure 

15.  Past  24-hr  change 
700-rab  height 

16.  Past  24-hr  change 
trough  height 

17.  Past  24-hr  change 
ridge  height 

18.  Past  48-hr  change 
element  #2 

19.  Past  48-hr  change 
element  #3 

20.  Past  48-hr  change 
element  #5 

21.  Past  48-hr  change 
element  //6 

22.  Past  48-hr  change 
element  #7 

23.  Past  48-hr  change 
(maximum  wind) 

24.  Past  48-hr  change 
pressure 

25.  Past  48-hr  change 
700-mb  height 

26.  Past  48-hr  change 
trough  height 

27.  Past  48-hr  change 
ridge  height 


in  derived 


in  derived 


in  derived 


in  intensity 


in  sea-level 


in  minimum 


in  700-mb 


in  700-mb 


in  derived 


in  derived 


in  derived 


in  derived 


in  derived 


in  intensity 


in  sea-level 


in  minimum 


in  700-mb 


in  700-mb 


+   2  deg  long 

+  10  m/deg  long 

+  10  m 

+  10  kt 

+   5  mb 

+  10  m 

+  10  m 

+  10  m 

+   2  deg  lat 

+  10  m/deg  lat 

+   2  deg  long 

+  10  m/deg  long 

+  10  m 

+  10  kt 

+   5  mb 

+  10  m 

+  10  m 

+  10  m 
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APPENDIX  F 

Indicated  Forecast  Potentials  of  the 
Derived  Analog  Elements 


Indicated  forecast  potentials  of  the  derived  analog 
elements  for  each  forecast  period  are  listed  below.   The 
indicated  potentials  are  expressed  as  a  percent  of  the 
total  number  of  possible  occurrences  in  which  the  trend 
of  decreasing  degree  of  matching  with  increasing  forecast 
error  occurred  for  each  element.   Refer  to  Appendix  E 
for  element  definitions. 


Analog  Element 
Numb  er 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 


Forecast  Potential  for 
Given  Forecast  Periods  (%) 


24-hr 

27 
20 
23 
20 
30 
40 
27 
29 
39 
21 
18 
29 
46 
32 
18 
29 
36 
32 
36 
28 
20 
36 
20 
28 
8 
24 
24 


48-hr 

31 
34 
34 
24 
31 
3  4 
21 
26 
41 
33 
26 
22 
37 
22 
37 
48 
33 
17 
25 
29 
17 
25 
21 
33 
4 
29 
29 


72-hr 

16 
20 
28 
16 
20 
23 
24 
25 
38 
21 
17 
25 
21 
21 
13 
21 
21 
14 
23 
18 

9 
27 
18 
18 

0 
23 

9 
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APPENDIX  G 
Selected  Predictors  for  each  Forecast  Period 

The  original  and  derived  analog  elements  that  were 
selected  as  predictors  for  each  forecast  period  are  listed 
below  along  with  their  associated  weight  factors. 

24-hr  Forecast  Period 


Predictor 


Weight  Factor 


Past  12-hr  change  in  sea-level  pressure 

Past  24-hr  change  in  intensity  (maximum  wind) 

Past  12-hr  speed  of  movement 

700-mb  trough  height-minimum  700-rob  height 
7  00-mb  trough  long-cyclone  long 

700-mb  ridge  height 
Past  24-hr  change  in 

700-mb  ridge  height-minimum  700-mb  height 

700-mb  ridge  lat-cyclone  lat 


7  00-mb  ridge  height 
thickness 


700-mb  ridge  lat 
Past  24-hr  change  in 
Past  48-hr  change  in 
Past  48-hr  change  in 

700-mb  ridge  heigh t -minimum  700-mb  height 

700-mb  ridge  lat-cyclone  lat 

700-mb  trough  long 
minimum  700-mb  height 
size 


48-hr  Forecast  Period 


Predictor 


.42 
.39 

.24 

.21 
.21 


1.1! 


.12 
.09 
.09 


1.09 


.00 
.00 
.00 


Past  24-hr  change  in  trough  height 
Past  24-hr  change  in 

700-mb  ridge  height-minimum  700-mb  height 

700-mb  ridge  lat-cyclone  lat 

700-mb  trough  height 

700-mb  ridge  height 

current  intensity  (maximum  wind) 

Past  24-hr  change  in  intensity  (maximum  wind) 

Past  24-hr  change  in  minimum  700-mb  height 

700-mb  ridge  lat 


Weight  Factor 
1.60 


1.37 


1 
1 
1 
1 
1 
1 


37 

27 
27 
23 
23 
13 
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4 8 - hr  Fo recast  Period  (continued) 


Predictor 

700-mb  ridge  height -minimum  700-mb  height 
700-mb  ridge  lat-cyclone  lat 

700-mb  trough  heigh t -minimum  700-mb  height 
700-mb  trough  long-cyclone  long 

Past  24-hr  change  in  700-mb  trough 
long-cyclone  long 

Past  24-hr  change  in  700-mb  ridge  height 
Past  48-hr  change  in  sea-level  pressure 
Past  12-hr  change  in  speed  of  movement 
Past  24-hr  change  in  direction  of  movement 

72-hr  Forecast  Period 


Weight  Factor 
1.13 

1.13 


1 
1 
1 
1 

1 


10 
10 
10 
07 
00 


Predictor 

Past  24-hr  change  in 
700-mb  ridge  height-minimum  700-mb  height 
700-mb  ridge  lat-cyclone  lat 

700-mb  trough  height 

700-mb  ridge  lat 

Size 

Past  12-hr  change  in  size 

700-mb  ridge  height-minimum  700-mb  height 
700-mb  ridge  lat-cyclone  lat 

Past  48-hr  change  in  thickness 

Past  24-hr  change  in  thickness 

Past  12-hr  change  in  sea-level  pressure 

Past  24-hr  speed  of  movement 


Weight  Factor 


1.81 

1.71 
1.52 
1.52 
1.33 

1.33 

1.29 
1.19 
1.14 
1.00 
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APPENDIX  H 
Independent  Test  Cyclones 

Independent  test  cyclones.   The  following  tropical 
cyclone  reports  were  selected  as  independent  test  cases 
for  the  analog  forecast  program. 


Name 


Cora 

Cora 

Polly 

Trix 

Trix 

Trix 

Bess 

Bess 

Wendy 

Wendy 

Wendy 

Hope 

Hope 

Marge 

Marge 

Opal 

Opal 

Opal 

Rita 

Rita 

Rita 

Alice 

Alice 

Alice 

Alice 

Alice 

Tess 

Tess 

Tess 

Cora 

Cora 

Cora 

Cora 

Cora 

Kim 

Kim 

Kim 


MO-DA-YR 


8-18-69 
8-20-69 
8-04-68 
8-24-68 
8-25-68 
8-26-68 
9-01-68 
9-03-68 
9-02-68 
9-04-68 
9-06-68 
8-05-67 
8-06-67 
8-24-67 
8-25-67 
9-05-67 
9-06-67 
9-07-67 
8-02-66 
8-03-66 
8-06-66 
8-27-66 
8-29-66 
8-31-66 
9-01-66 
9-02-66 
8-12-66 
8-13-66 
8-14-66 
8-31-66 
9-01-66 
9-02-66 
9-04-66 
9-05-66 
8-04-65 
8-05-65 
8-06-65 


Time 

Initial 

(GMT) 

Intensity 

(kt) 

1200 

65 

0600 

80 

0600 

55 

0600 

35 

1200 

40 

0600 

50 

0000 

40 

0600 

50 

0000 

105 

0000 

85 

0600 

60 

0600 

45 

0000 

55 

0600 

30 

0000 

40 

0600 

110 

0600 

115 

0600 

P5 

1200 

25 

1200 

35 

0600 

50 

0000 

55 

1800 

75 

1800 

100 

1200 

130 

1200 

130 

1200 

25 

1200 

35 

1800 

55 

0600 

50 

1200 

100 

1800 

130 

0000 

130 

0600 

140 

0600 

25 

0600 

50 

0600 

60 
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Name 


Lucy 

Lucy 

Lucy 

Lucy 

Lucy 

Olive 

Olive 

Olive 

Kathy 

Kathy 

Kathy 

Kathy 

Kathy 

Kathy 

Kathy 

Kathy 

Kathy 

Marie 

Marie 

Marie 

Sally 

Sally 

Sally 

Sally 


MO-DA-YR 


8-16-65 
8-17-65 
8-19-65 
8-20-65 
8-21-65 
8-29-65 
8-30-65 
8-30-65 
8-12-64 
8-13-64 
8-14-64 
8-15-64 
8-17-64 
8-18-64 
8-19-64 
8-20-64 
8-22-64 
8-15-64 
8-15-64 
8-16-64 
9-04-64 
9-05-64 
9-06-64 
9-07-64 


Time 

Initial 

(GMT) 

Intens  ity 

(kt) 

0000 

85 

0000 

120 

0000 

120 

0600 

150 

0600 

120 

0000 

150 

0000 

130 

1800 

110 

0600 

35 

0600 

80 

1200 

80 

1800 

60 

0000 

75 

0600 

110 

1200 

110 

1800 

100 

0000 

85 

0600 

30 

1800 

30 

0600 

30 

0000 

55 

0000 

85 

0600 

130 

1800 

160 
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